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We propose a method to measure the normal fraction of a two-dimensional Bose gas, a quantity 
that generally differs from the non-condensed fraction. The idea is based on applying a spatially 
oscillating artificial gauge field to the atoms. The response of the atoms to the gauge field can 
be read out either mechanically from the deposited energy into the cloud, or optically from the 
macroscopic optical properties of the atomic gas. The local nature of the proposed scheme allows 
one to reconstruct the spatial profile of the superfluid component; furthermore, the proposed method 
does not require having established thermal equilibrium in the gas in the presence of the gauge 
field. The theoretical description of the system is based on a generalization of the Dum-Olshanii 
theory of artificial gauge fields to the interacting many-body context. The efficiency of the proposed 
measurement scheme is assessed by means of classical field numerical simulations. An explicit atomic 
level scheme minimizing disturbing effects such as spontaneous emission and light-shifts is proposed 
for 87 Rb atoms. 
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I. INTRODUCTION AND SUMMARY OF MAIN 
RESULTS 

One of the most striking features of degenerate Bose 
gases in two dimensions is the possibility of having a 
superfluid behavior in the absence of a macroscopically 
populated Bose-Einstein condensate. The transition 
to the superfluid state is of the Bcrezinskii-Kostcrlitz- 
Thouless (BKT) type, characterized by a sudden jump of 
the superfluid density from to the universal value 4 (in 
units of the inverse square of the dc Broglic thermal wave- 
length) , independent from the details of the system [J, Q ■ 
At the transition point, the asymptotic behavior of the 
field correlation function changes from an exponential to 
a power-law decay at large distances. In contrast to the 
three-dimensional case, superfluidity is then not related 
to the appearance of a macroscopically occupied Bosc- 
Einstcin condensate in the thermodynamic limit. 

Pioneering experiments have addressed the mechani- 
cal properties of two-dimensional layers of liquid Helium 
adsorbed on a substrate [|[ and have characterized the 
universal jump of the superfluid fraction at the BKT crit- 
ical point. On the other hand, liquid Helium experiments 
have limited access to the momentum distribution and 
the correlation functions of the fluid. The situation of 
ultracold atom experiments is almost the opposite: ev- 
idence of the BKT transition has been obtained from 
the coherence functions [|J , the number of observed vor- 
tices [f|, and the density profile after time-of- flight @, 
while the macroscopic mechanical properties of the fluid 
have not been characterized yet. 

Quite some effort has been recently devoted to the con- 
ceptual problem of how to experimentally detect genuine 
superfluidity in a quantum gas of ultracold atoms and 
not simply Bose-Einstein condensation Q. A possibility 
explored in Q is to look at the response of a gas in a 



toroidal trap to a static azimuthal artificial gauge field: 
a spectroscopic signature is proposed which should pro- 
vide direct information on the total superfluid mass of the 
system. A different strategy proposed in Q consists of 
looking at the evolution of the equilibrium density profile 
of a trapped gas when it is set into rotation. 

In the present paper we propose two experimental pro- 
tocols to measure the normal fraction of a gas in a local 
way, so to extract its spatial dependence in a trapped 
geometry. This feature is most relevant for atomic sam- 
ples, as the superfluid core coexists with an external ring 
of normal gas In particular, the proposed diagnostic 
technique does not require to relate experimental obser- 
vations after time of flight to in-trap quantities. Fur- 
thermore, in contrast to @ our technique does not re- 
quire thermodynamic equilibrium in the gas in presence 
of rotation pd| and may be applied to more general, non- 
equilibrium conditions. 

The basic idea of our proposal is based on the definition 
of normal and superfluid fractions of a quantum fluid in 
terms of its current response to a transverse gauge field 
in the low- frequency and long- wavelength lim it [fa . [l6| . 
A spatially oscillating artificial gauge field p7H20| with a 
spatially localized envelope can be applied to the atomic 
gas using a suitable combination of laser beams. The 
response of the fluid to the gauge field can be detected 
either mechanically or optically. In the former case, one 
has to measure the amount of energy that is deposited in 
the atomic gas at the end of a suitable temporal sequence 
of gauge field. In the latter case, one can observe e.g. the 
phase shift that is experienced by the laser fields while 
crossing the atomic cloud. 

The structure and the main results of the paper can 
be summarized as follows. In Sec|H]we review the defini- 
tion of the normal and superfluid fractions that wc adopt 
throughout the whole paper. A strategy to generate the 
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artificial gauge field with the suitable spatial geometry 
is presented in Sec lIIII using three laser beams, namely a 
coupling beam and two probe ones. 

The first method to measure the normal fraction is by 
mechanical means. It is discussed in Sec lIVt a pulse of 
spatially modulated gauge field is suddenly applied to 
the gas and then slowly switched-off according to an ex- 
ponential law in time. An analytical calculation within 
the linear response theory and local density approxima- 
tion shows that the energy that is deposited in the gas 
at the end of the gauge field pulse is indeed proportional 
to the normal (total) density in the small fraction of the 
gas where the probe beam is focused, if the spatial mod- 
ulation of the gauge field is orthogonal (parallel) to the 
coupling beam direction. Some of the complications that 
naturally occur in experiments are then numerically in- 
vestigated: a Bogoliubov theory is used to assess the con- 
ditions to be imposed to the geometry and the temporal 
duration of the pulse; a classical field model is used to as- 
sess the conditions on the amplitude of the gauge field for 
the linear response theory to be valid. A main difficulty 
appears to be the relatively small amount of energy that 
can be deposited in the gas before nonlinear couplings 
become difficult to extrapolate out: the resulting figure 
is on the order of 1% of the total energy of the gas, which 
is however not far from the sensitivity of state-of-the-art 
thermodynamic measurements of the energy (2l| - [25j . 

The second method to measure the normal fraction is 
by optical means. It is discussed in Sec IV Al The spa- 
tially modulated gauge field is imposed using the same 
laser beam configuration in a continuous wave regime to 
generate a stationary current pattern in the gas. This 
current pattern can be read out from the phase shift ac- 
cumulated by the same probe laser beams after crossing 
the atomic cloud. Analytical calculations show that this 
phase shift is indeed proportional to the normal (total) 
density if the spatial modulation of the gauge field is or- 
thogonal (parallel) to the coupling beam direction. For 
realistic configurations, the amount of the phase shift is 
anticipated to be of the order of a fraction of 10 -3 , i.e. 
small but still appreciable with present-day optical tech- 
niques. 

Another optical set-up that is able to provide quanti- 
tative information on the normal fraction is discussed in 
Scc lVBt analytical calculations are used to relate the an- 
gular distribution of the scattered light off a single probe 
beam to the normal fraction of the gas. The main advan- 
tage of this last configuration is the rapidity of the mea- 
surement but the limiting factor is the relatively small 
number of photons (of the order of a few units) that are 
expected to be scattered in the useful directions. 

In Appendix [X] we discuss how to cope with all those 
spurious effects that arise from a realistic configuration 
of atomic levels and laser fields, in particular sponta- 
neous emission and the mechanical effect of undesired 
light-shifts; even though the underlying concepts are gen- 
eral, the discussion is mostly focused on the most promis- 
ing case of 87 Rb atoms. The generalization of the Dum- 



Olshanii theory of artificial gauge fields to the many-body 
context is reported in Appendix [Bj this development is 
required to put the gauge field concepts on firm ground 
in presence of atomic interactions and to evaluate in a 
rigorous way the optical response of the atoms to the 
combined coupling and probe beams. Appendix [Q gives 
more details on the analytical derivation of the deposited 
energy and clarifies some issues related to the local den- 
sity approximation. The framework for calculating the 
deposited energy using the Bogoliubov theory within the 
linear response regime is discussed in Appendix [D] The 
last Appendix [E] discusses issues related to the statisti- 
cal noise on the deposited energy both in the numerical 
calculation and in an actual experiment. In particular, 
it shows how useful information on the normal fraction 
could be extracted from the noise if a sufficiently pre- 
cise determination of the initial energy was possible. We 
conclude in section |VT1 



II. 



DEFINITION OF SUPERFLUID AND 
NORMAL FRACTIONS 



Our proposal to quantitatively assess the superfluidity 
of the two-dimensional atomic gas is based on the tra- 
ditional definition of the normal fraction /„ in terms of 
the response to a transverse gauge field coupling to the 
atomic current operator (lH [la ] . The Hamiltonian giving 
the coupling of the matter current to an arbitrary vector 
potential A(r) is 
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with the current operator defined as usual as 
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in terms of the bosonic field operator (j> for the two- 
dimensional gas. For a spatially homogeneous system, 
the linear response susceptibility relating the average cur- 
rent [27[ to the applied gauge field can be easily written 
in momentum space and frequency domain as 



(j)(q,w) = x(q,w) A(q,w) 



(3) 



If the system is also invariant under reflection with 
respect to the direction of q, the susceptibility ten- 
sor x(q, ui) turns out to be diagonal in the longitudi- 
nal/transverse basis with respect to q, with diagonal ma- 
trix elements Xl,t{h, &>), respectively 

For a system of surface density p, the normal fraction 
f n of the system is then defined as the low-momentum, 
low-frequency limit of the susceptibility to transverse 
gauge fields: 



/„ = lim lim — xr(q, w). 

<2->0cj-s-0 p 



(4) 
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Note that the order of the limits is here important. A 
well-known sum rule based on gauge invariance imposes 
that the same limit for the longitudinal susceptibility 
Xi(q, oj) gives exactly unity, 

JJJ 

1 = lim lim — XL(q,w). (5) 

q->0 Ul— >0 p 

The definition ((4]) can then be extended to large but 
finite systems using the standard local density approxi- 
mation. 



III. HOW TO GENERATE THE GAUGE FIELD 

We consider a three-dimensional gas of bosonic atoms 
in a strongly anisotropic, pancake-shaped trap. The axial 
confinement frequency co z is much higher than the one ui\\ 
along the xy plane; both the temperature T (times the 
Boltzmann constant ks) and the chemical potential /i 
of the gas are assumed to be smaller than Hui z . In this 
regime, the gas will be eventually described in terms of 
a two-dimensional Hamiltonian. 

Building on an idea originally introduced in [l7j], an 
artificial gauge field coupling to the atomic current can 
be obtained by illuminating the atoms with several laser 
beams with suitably chosen frequencies, wavevectors, 
and waist profiles. Several schemes to generate artifi- 
cial gauge fields for neutral atoms have been proposed in 
the last years [l7l - [l9| . The last proposal [l9| was recently 
implemented on an atomic Bosc-Einstein condensate: for 
sufficiently strong gauge fields, a, disordered ensemble of 
vortices appeared in the gas [20| . 
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FIG. 1: Scheme of the set-up under consideration. Left 
panel: generic sketch of the A configuration of atomic lev- 
els and laser beams involved in the optical processes. Center 
panel: view from above of the two-dimensional atomic pan- 
cake (lying within the dashed circle) and of the laser beams 
(the hatched disk is the spot of the probe beams in the xy 
plane). Right panel: side view. k c and are the wavevec- 
tors of the coupling and probe beams, respectively. In prac- 
tice, kp = (kp — q 2 /4) 1,/2 e z ± (q/2) e x with q <C k p , where 
is the unit vector along axis i. 

In the present paper, we shall focus our attention on 
the level configuration shown in FigfTJ Three internal 
atomic levels in a A configuration are connected by three 
laser fields according to the sketch given in the left panel 



of FigfT] a coupling beam resonantly drives the \b) —> |e) 
atomic transition, while a pair of probe beams resonantly 
drive the |o) —> |e) transition with essentially the same 
detuning. The artificial gauge field originates from the 
spatial and temporal dependence of the resulting opti- 
cally dark state [13] • All other atomic states are assumed 
to be far-off resonance. A discussion of their effect in the 
specific case of 87 Rb atoms is given in the Appendix [A] 
A summary of the suggested experimental parameters for 
this specific atomic species is given in Table HI 

The geometrical arrangement of the laser beams is 
sketched in the central and right panels of FigfT] The 
continuous- wave control beam propagates along the y di- 
rection with a wavevector k c and a carrier frequency uj c 
close to resonance with the \b) — > |e) transition, with a 
detuning d, and has a peak Rabi frequency Its waist 
profile is much wider than the size of the atomic cloud, 
so that it can be safely approximated by a plane wave. 

The two probe beams share the same carrier frequency 
Up close to resonance with the \a) — > \e) transition. The 
carrier frequencies of the coupling and probe beams are 
chosen exactly on resonance with the Raman transition 
| a) — > |e) — > \b), that is uj p — lu c = lui, — ui a . The probe 
beams impinge on the atomic cloud with wavevectors 
k^ ~ k p e z ± q/2 close to the z direction and sym- 
metrically located with respect to it. The difference 
q = k+ — k~ lies along the xy plane and is in magni- 
tude q <C k v — ujp/c. The probe beams spatial profile is 
taken to be a Gaussian with a waist w, centered at ro in 
the z = plane. 

In what follows, we will need the condition w 3> <z _1 , 
where q is small enough for the limit in Eq.Q to be 
almost reached. In the numerical examples to come, this 
requires q£ < 1, where £ is the healing length of the 
gas. This automatically shows that the concept of normal 
fraction can be used only when the portion of the gas 
within the waist is "macroscopic" [Hj], here w ^> £, which 
was physically expected. At the same time, the waist w 
is assumed to be much smaller than the atomic cloud 
radius, so to allow for a local measurement of the normal 
fraction. The peak Rabi frequencies of the two probe 
beams are (t) , respectively (see Appendix [B] for the 
precise definition of the Rabi frequencies). The spatial 
dependence of the Rabi frequencies of both the coupling 
and the probe beams is then summarized by the following 
expressions [29l |. 

n c (r,t) = ^e-^'e*" (6) 
xe -l( x - x o) 2 +(y-yo) 2 ]/w 2 ^ ^ 

where we have allowed for the coupling beam to have a 
small detuning A c from Raman resonance on top of its 
carrier frequency at tu c . On the contrary, H p does not 
have a time dependent phase factor, but only contains 
the square root of a purely real non-negative switch-on 
and switch-off function f(t), common to the two probe 
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beams: 



il p (r,t) = il° p (r)[f(t)] 1 ^. 



(8) 



In what follows, we shall restrict our attention to the 
non-saturating regime |fi c |, |f2p | <C |<5 + ir/2|, where T is 
the decay rate of |e) due to spontaneous emission, and 5 
is the common detuning of the probe and control beam 
carrier frequencies from the transitions |o) — > |e) and 
| b) — > |e) respectively. We shall also concentrate on the 
limit | <C |O c | where the structure of the gauge field 
is the simplest. 

As the transitions driven by the probe and control 
beams share the excited state |e), there exists for each 
spatio-temporal coordinates (r, t) an internal non-coupled 
state, for which the two excitation channels interfere de- 
structively. In terms of the local Rabi frequencies fi p (r, t) 
and f2 c (r, t), this non-coupled state reads 



\NC{T,t)) 



\a)-n p (r,t)/n c {r,t)\b) 
(l + |n p (r,t)P/|n c (r,t)|2)V2- 



(9) 



Adiabatically eliminating the excited state |e), one sees 
that the bright orthogonal state, the so-called coupled 
state, 



\C(v,t)) 



[n p (r,t)/n c (r,t)]*\a) + \b) 

(i + \n p (r,t)\y\n c (r,t)\ 2 y/* 



(10) 



is separated from \NC(r,t)) by a (complex) energy gap 

h[S (r,t)-.r (v,t)/2] = 4{S + ir/2) > ( n ) 

where 6' and V are the light-shift and the decay rate of 
the coupled state. If the energy gap is large enough as 
compared to both the motional coupling between \NC) 
and \C) due to the spatio-temporal dependence of f2 c and 
f2^ [30| , and to the quantum of oscillation tko z along the 
tightly confined z direction, we can restrict the dynamics 
to the \NC) internal state. 



Generalizing the single-particle theory of [17[ to the 
many-body context in Heisenberg picture, one gets to an 
effective Hamiltonian for the component <p3u(r,t) of the 
three-dimensional atomic field operator in the (spatially 
and temporally-dependent) non-coupled state \NC(r, t)), 

$ 3 D(v,t) = (NC(r,t)\a) $ a (r,f) + (NC(r,t)\b) $ 6 (r,t) 

(12) 

in the simple form [31| : 
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U 3D {r) + W 3D (r,t) 



9.315 



-j 3Z5 (r) ■ A 3D (r,t) + ^g{r,t)4>l D (f)l D <f>3D^3D> (13) 

where U 3 d is the three-dimensional trapping potential 
(supposed to be common to the internal states \a) and 
| b)) and the vector gauge potential 



A 3D (r,t) 
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couples to the three-dimensional atomic current operator 
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J3i5(r) 
and the scalar potential 



^3i5( r ) V ^3i5(r) - h.c. 



(15) 



W 3D (r,t) = - l -[(NC(r,t)\ [d t \NC(r,t))} - c. 

+ ¥- E idrANC(r,t)\}[d n \NC(r,t))} (16) 
2m . 



couples to the three-dimensional density 
^3/5 (r) = 4>l D (r) 4>sd(t). 



(17) 



The derivation of the Hamiltonian (jTB")) is based on 
the Quantum Stochastic Differential Equations formal- 
ism [32( ; the details are given in the Appendix [B] 

The spatial and temporal dependence of the weights 
of the non-coupled state © in the internal states \a) 
and | b) reflects into a similar variation of the coupling 
constant describing the atomic interaction within the in- 
ternal state | AC): 



9ZD{v,t) 



\tt c \ 4 g aa + 2\n p {v,t)\ 2 \n c \ 2 9ab + \n p \± gbb 
(|fi c P + |f2 P M)P)2 



(18) 

where the coupling constants g aa , g ab and g bb originate 
from the a — a, a — b and b — b elastic s-wave interac- 
tions. In what follows, we shall be interested in isolating 
the response of the system to the gauge field A 3 d- To 
this purpose, it will be useful to minimize the effect of 
all unwanted couplings to the density introduced by the 
scalar potential W 3 d and by the spatio-temporal depen- 
dence of the interaction constant g 3 o- This latter effect is 
minimized if one chooses states a, b with similar scatter- 
ing properties g aa ~ g ab ~ g bb . In the limit |f2 p /f2 c j < 1, 
one simply needs to have g aa ~ g ab , as it is assumed from 
now on. 

The atomic field along z is assumed to be frozen in the 
ground state of the harmonic confinement of wavefunc- 
tion (j)o(z). This allows to express the three-dimensional 
bosonic field (f> 3 jj in terms of the bosonic field for a two- 
dimensional gas, setting 



cf) 3D (x,y,z) = <f) Q (z)(j>(x,y). 



(19) 



Correspondingly, the two-dimensional coupling constant 
g has the expression 



9 = —9 = 
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(20) 



ho 



in terms of the three-dimensional coupling constant 
93D = 4:irh 2 a 3 D/m, where a 3 o is the s-wave scattering 
length in the state \a), and the size a{; o = \Jh/mw z of 
the ground state along z. Note that the dimensional re- 
duction (TT9")) does not require to be in the Lamb-Dickc 
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limit fc PjC a^ <C 1, it is sufficient that the typical energy 
per particle corresponding to the various terms in % is 

< hio z . 

The effective two-dimensional gauge and scalar poten- 
tials A and W then result from an average of the Hamil- 
tonian (fT5|) over the motional ground state along z. In- 
cluding in (p~4)) and (fT6|) the explicit form of the beam 
profiles and restricting ourselves to zcroth order in the 
small parameters g/fc PiC , l/wk PtC , and to second order 
in $7^ /fi c , the resulting two-dimensional gauge potential 
turns out to be directed along the y axis and to have the 
form 



A y (r) ~ Hk c 
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\n+\ 2 + \n-\ 2 + [n+n-*e^-' + c.c.] 
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-2|r-r | 2 /tt> 2 _ ^21) 



To the same level of approximation, the scalar potential 
has the form 



W(r) 



2m 



HA, 



(22) 



which can be made to vanish by choosing a detuning A c 
that exactly compensates the recoil of the atoms after 
the Raman process a — > e — > b: 



W = for A r = - 



ml + k 2 p ) 



2m 



(23) 



It is worth pointing out that the temporal derivative 
df(t)/dt of the real switch-on and switch-off function f(t) 
of the probe beam Eq. ((SJ) has an exactly vanishing con- 
tribution to the d t term in the right-hand side of Eq. (fTBf , 
so it does not contribute to the scalar potential W to all 
orders in Q p /Q c [33|. The same conclusion holds for the 
gauge field, since there is no temporal derivative in (|T4"|) . 

After expansion of the squared modulus as done in the 
second line of (f2~Tj). two kinds of terms are immediately 
identified: (i) a slowly varying Gaussian term of size w 
and peak amplitude | 2 + |fi~ | 2 that follows the laser 
envelopes, and (ii) an oscillating term at wavevector q 
with a Gaussian envelope of size w and peak amplitude 
|fi+Q~|. This spatially modulated term is indeed the one 
that we need to probe the normal fraction of the gas ac- 
cording to the definition (UJ: when q is taken along the 
x axis (y axis), it provides an almost purely transverse 
(longitudinal) contribution to the gauge field A. On the 
other hand, the slowly varying term always includes both 
longitudinal and transverse vector field components. Ex- 
perimental procedures to subtract the effect of this un- 
wanted term will be discussed in the next sections. 



IV. DEPOSITED ENERGY MEASUREMENT 
A. General idea 

In this section we shall present a method to extract 
the value of the normal fraction from a measurement of 
the energy that is deposited in the system by a suitably 
designed gauge field sequence. The coupling beam is as- 
sumed to be always switched on. On the other hand, both 
probe beam intensities |f2p | 2 are varied in time according 
to the (dimcnsionless) real envelope function /(t). This 
function is chosen to be for t < and to rapidly tend 
back to at long times. 

As already mentioned, we assume that the atomic in- 
teraction constants satisfy g aa ~ g a b [HI- As soon as 
qw 3> 1, the deposited energy is the sum of two inde- 
pendent contributions AEi^ corresponding to the de- 
composition (|2~T|) of the gauge field as the sum of a non- 
modulated term and a modulated one at the wavevector 
q. Using standard linear response theory within the lo- 
cal density approximation as discussed in the Appendix 
[Cl the contribution AE2 of the modulated term can be 
written in the simplified form 



t 2 ( fgaugc \ 2 f + °°du; 

4"' 



AE 2 ~ -w 2 ' ga 



f — Lj\f(Lj)\ 2 lm[xyy{q,uj)} 

J-oo * 



(24) 

where we have introduced the amplitude of the spatially 
modulated part of the gauge field, 



^gauge — 2 Hk c 



|n+n-|(t = o+) 



(25) 



A similar expression for the contribution AE\ of the non- 
modulated term is given in the Appendix [C] as (|C8|) . The 
susceptibility Xyy appearing in (|24[) is the one (as defined 
in section [TTJ) of the fictitious spatially homogeneous two- 
dimensional Bose gas that approximates the state of the 
trapped gas around the center ro of the probe beams' 
spot. 

The expression ([24| for the deposited energy involves 
the imaginary part of that susceptibility, while the nor- 
mal fraction §4§ involves the real part. To relate the two, 
one can make use of the well-known Kramers-Kronig re- 
lation of linear response theory, 



lim Rc[x(q,w)] 

uj— >0 



+00 



duj' Im[x(q,o/)] 

7T U>' 



For a suitably chosen envelope of the form 



/(*) 



0(t), 



(26) 



(27) 



where Q(t) is the Heaviside step function, the Fourier 
transform is equal to f(u>) = i/(u) -\- 17), and the integral 
in (f2"4"]) indeed reduces to the real part of the susceptibility 
([26]) in the 7^0 limit. 

As a consequence, the deposited energy AE 2 for small 
q perpendicular (parallel) to k c can be related to the nor- 
mal (total) density p n (p) of the trapped gas at position 
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AEo ~ — 



7T W 
4 
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Pn(r ) 



for q 1 k c 
for q || k c 



(28) 



For intermediate angles a between q and k c , AE2 is pro- 
portional to p cos 2 a + p, 



, sin 2 a 



In an actual experiment, the undesired contribution 
AEi can be eliminated by noting its independence on 
the relative orientation of q and k c [see Eqs. (|C8[) and 
(|C9j) ]. as well as its different dependence on the probe 
amplitudes ftp, proportional to [|f2+ 1 2 + |fi~| 2 ] 2 rather 
than |f2+f2~| 2 . By measuring the deposited energy for at 
least two different values of the f2+ ratio, one is able 
to isolate the relevant contribution 
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B. How fast is the q, 7 — > limit reached ? 

An important point in view of experiments is to char- 
acterize how small the wavenumber q and the switch- 
off rate 7 of the gauge field spatial modulation have ac- 
tually to be taken to obtain a quantitatively accurate 
measurement of the normal fraction /„ . To answer this 
question, we consider in this subsection the simplest case 
of a spatially homogeneous two-dimensional system in a 
square box of size L with periodic boundary conditions, 
excited by the gauge field in a plane wave form, i.e. in 
the limit w — > 00. We also limit ourselves to the case of 
a transverse gauge field with q = qe x perpendicular to 
k~ k r G11 , 
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(r,t) = e-^e(t)e y ^(e 
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(29) 



The deposited energy at the end of the gauge field se- 
quence can be evaluated by means of the Bogoliubov 
theory of dilute Bosc gases. The main steps of the calcu- 
lation are sketched in the Appendix ID] The final result 
reads 



;c ^ 2 _ fefi 

_ J n 



(30) 



in terms of the wavevector- and 7-dependent effective 
normal fraction /° ff : 



f cff = — \ ~" 



h 2 kl 



N f — ' m 

k^O.-q 



^ [ - " k+q fc „ (C/kC/k+q - V k Fk+q) 2 



+ 



L ek+q - ek - ifry 

1 + «k + ™k+q 



ek - ifvy 



(C/kVk+q - ^kC/k+q)^ 



(31) 



Here N is the total particle number, p = N/L 2 is the 
surface density, 
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FIG. 2: Quantum Bogoliubov prediction (|31|) for the effec- 
tive normal fraction /® ff . In (a), dependence of /® ff on the 
switch-off rate 7 of the gauge field, for three different values of 
the gauge field spatial modulation wavevector q = (2n/L) e x 
(solid, q£ ~ 0.1), q = 7 x (2Tv/L)e x (dotted, g£ ~ 0.7) 
and q = 10 x (2tt/L)b x (dashed, q£ ~ 1). System pa- 
rameters: square box (with periodic boundary conditions) of 
size L/f ~ 63 containing N ~ 40000 particles with interac- 
tion constant g — 0.1h 2 /m at a temperature T/Ta = 0.1. 
ksTd = 2nh 2 p/m is the degeneracy temperature. The heal- 
ing length £ of the gas is defined by ti 2 /(m£ 2 ) = pg and the 
Bogoliubov sound velocity is defined as c„ = (pg/m) 1 ^ 2 . The 
dashed curve corresponds to the same value of g£ ~ 1 as 
in Fig(4] and the vertical dotted lines indicate the values of 
~//(c s q) considered in that figure. The thin horizontal line 
is the prediction (|34[) of the quantum Bogoliubov theory in 
the thermodynamic limit. Interestingly, for the considered 
atom number, the finite-size thermodynamic formula Eq. (|39[) 
evaluated with the quantum Bogoliubov theory gives already 
the same value at the scale of the figure [37j|. The quadratic 
rather than linear dependence of f„ s on 7 for small values of 
7 is a finite-size effect. The fact that in (a), the values of f„ a 
for the dotted curve are close to the arithmetic mean of the 
solid and dashed curves, for common values of ^/(c a q), sug- 
gests that fn B is roughly an affine function of (q£) 2 for fixed 
7/(c s g). This guess is substantiated in panel (b) where the 
dependence of the effective normal fraction on the modulation 
wavevector q = qe x is shown for y/(c s q) — 0.2 [the crosses 
are the Bogoliubov predictions and the line is a linear fit of 
/^ ff as a function of q 2 ]. 



7 



is the usual Bogoliubov dispersion relation and the am- 
plitudes of the Bogoliubov modes satisfy 



Uu + 14 = 



1 



h 2 k 2 /2m 
H 2 k 2 /2m + 2pg 



1/4 



(33) 



The nk are the thermal mean occupation numbers of Bo- 
goliubov modes, rtk = l/{exp(ek/kBT) — 1]. 

The thermodynamic limit L — > oo at fixed density 
N/L 2 can be worked out analytically by first taking the 
7 — > limit and then the q — > limit in the expres- 
sion (|31|) for /° ff . In this way one recovers the usual 
Bogoliubov expression for the normal fraction, which in 
dimension two reads: 



i r d 2 k h 2 k 2 

fn = -pj (2^~ ( ~ 9 - nk) - 



(34) 



For a finite size system, the dependence of /^ ff on 7 and q 
is explored in FigJ5] For the smallest non-zero wavevector 
value allowed by the chosen quantization box, the relative 
error on f° s is already as small as 10% for -f/c s q = 0.15, 
where c s = {pg/m) 1 / 2 is the Bogoliubov sound velocity. 
For 7/(0,5(7) fixed, in particular to the relevant value 0.2, 
and over the relevant range qt; < 1, we have also found 
that is to a good approximation an affinc function of 
(<?£) 2 > see Fig [2b, which was not obvious. 

Another interesting result of Bogoliubov theory ap- 
plied to our system is a sufficient condition on the am- 
plitude of the gauge field to be within the linear re- 
sponse regime. To this purpose, we can write the equa- 
tions of motion for the Bogoliubov mode operators 6k 
in the interaction picture in the presence of the time- 
dependent gauge field, and impose that the amplitude 
change is small as compared to the initial value. The 
expression of this amplitude change, that we shall not 
give here, involves as usual energy denominators such as 
£k - £k± q + ifhy. 

In the limit 7 — > 0, the condition of linear response is 
thus most stringent for modes of wavevector k such that 
£k = £k±q, where the real part of the energy denominator 
can vanish. For the maximal value of k set by the thermal 
occupation, this leads to the sufficient condition 



Cgaugo 



(mk B T d y/ 2 



< 



1/2 



2% 
kuT d 



(35) 



in terms of the degeneracy temperature ksTj = 
2irh 2 p/m. This naive argument is however not able to de- 
termine to which extent this condition is actually neces- 
sary. This would require a higher order calculation which 
falls beyond the scope of the present work [53] . 



C. Numerical investigation 

To further assess the validity and accuracy of our pro- 
posed scheme we have performed full scale numerical sim- 
ulations of the response of a two-dimensional Bose gas at 



finite temperature to the complete gauge field (f2Tj) , in- 
cluding the Gaussian envelope of the gauge field and a 
circular well trapping potential. A very useful tool to this 
purpose is the classical field model developed and applied 
in a number of recent works [lj], ■ For this model both 
the thermal equilibrium state and the temporal dynamics 
can in fact be easily addressed with standard numerical 
tools and provide reliable results for the physics of the 
degenerate Bose gas. 

We consider a classical (c-number) complex field de- 
fined on a square grid. The real space lattice constant b 
is chose n in terms of the thermal de Broglie wavelength 
A = v / 2nh 2 /mknT as b/X = A / 7r /( 4 C)- Th e specific 
value C — 0.80 [39( of the numerical coefficient is such 
that the classical field model correctly reproduces the to- 
tal number of non-condensed particles for an ideal gas 
at zero chemical potential in the thermodynamic limit. 
This choice corresponds to setting the ultra-violet mo- 
mentum cut-off fcmax = ir/b at h 2 k 2 alix /2ni = C,ksT in 
the classical field theory. 

In the canonical ensemble, the thermal probability dis- 
tribution for the interacting classical field follows a Boltz- 
mann law 5(||*|| 2 — N) cxp(— E[ip]/ksT) with the norm- 
squared 1 1* || 2 = 6 2 l*( r )l 2 fixed to the total atom 
number, and the discrete Gross-Pitaevskii energy func- 
tional [4f| given by: 



E[V] =b 2 J2^* 



-^A + U(r) + y - |*| 2 
2m 2 



*, (36) 



where U(r) is the trapping potential seen by the bidimen- 
sional gas. This probability distribution can be sampled 
by the long time limit of a Ito stochastic differential equa- 
tion with drift terms and a noise term fill [42j , including 
projectors in order to keep the norm-squared constant, 
\m\ 2 = N: 



ft2 AT 



2m 



U ^ + g\V\ 2 V 



Vk^T M-l 
— I — Q*a£ - 



2N 



k B Tdrm, (37) 



where M is the number of grid points, Qq, is the projector 
onto the subspace orthogonal to the classical field ty, d£ 
is a complex Gaussian, zero-mean, delta-correlated and 
temporally white noise such that dt;dt; = and 



dC(r l )d^(r J ) = drS r 



(38) 



Our numerical procedure simply consists in first gener- 
ating a number n roa i of independent wavefunctions dis- 
tributed according to the thermal Boltzmann law with 
energy functional (|36|) and then to let them evolve in 
real time according to the discrete time-dependent Gross- 
Pitaevskii equation including the gauge potential (|2"Tj) 

As a first application of the classical field method, we 
have determined the normal fraction of a spatially homo- 
geneous two-dimensional interacting gas at thermal equi- 
librium. This is done using the thermodynamic formula, 
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FIG. 3: Classical field simulation of the normal fraction /„ for 
a finite size two-dimensional, spatially homogeneous interact- 
ing Bose gas as a function of temperature. The calculation has 
been performed using the thermodynamic expression (|39|) for 
the normal fraction. The coupling constant is g = 0.1 h 2 /m. 
The different curves refer to simulations performed with dif- 
ferent system sizes L, resulting in different numbers of grid 
points M — 16 2 (green, dotted), 32 2 (red, dashed), 64 2 (black, 
solid) [the grid spacing is b ~ 0.99A, where A is the thermal 
de Broglie wavelength, see text] . The number of classical field 
realizations is n rca i = 1000. The dot-dashed line is the clas- 
sical field prediction for the normal fraction f n of an ideal 
gas in the thermodynamic limit: as discussed in the text, the 
decrease at high temperatures is an artifact of the classical 
field model. The blue circle indicates the result of a numeri- 
cal simulation of the deposited energy scheme as in FigUfor a 
homogeneous system with an infinite beam waist to = oo (in 
which case we simply dropped the non-modulated term in the 
gauge field), a number of grid points M = 32 2 , a gauge field 
modulation wavevector q = (2ir/L)e x ~ (0.2/£) e x , and in- 
cluding a careful extrapolation of e gaU gc — > (numerics down 
to e g augc/(77ifcsrd) 1/ ' 2 = 0.01) and of 7 — > (numerics down 
to 7/(c s g) = 0.05 where c s = (pg/m) 1 ^ 2 is the Bogoliubov 
sound velocity). 



applicable to finite size systems with periodic boundary 
conditions, 



/•the 



NmksT 



(39) 



which involves the thermal variance of the total momcn- 



in the case of an ideal gas, for which one can show that 
/„ = 1 - T/T d + 0(e- T <*/ T T/T d ) in the thermodynamic 
limit. 

The experimental estimation of the normal fraction 
obtained by the deposited energy method discussed in 
Scc lIV Al is simulated in FigHJ The value of the deposited 
energy is extracted from the classical field simulation by 
taking the energy difference at the end of two evolutions 
using the same value of |f2+ | 2 +|f2~ | 2 but different relative 
magnitudes of il^ , fl~ = and | f2+ 1 = | | respectively. 
This protocol aims at isolating the effect of the spatially 
modulated gauge potential: in the linear response limit, 
it is able to provide the exact value of AE2 alone. The 
effective normal fraction is then extracted from the de- 
posited energy via (|28|) , and it is plotted in Fig HI as a 
function of the gauge field amplitude, for different val- 
ues of the switch-off rate 7/(c s g) = 0.4, 0.2, 0.1 (in black, 
in red, in blue, from bottom to top) [c s — {pg/m) 1 / 2 is 
the Bogoliubov sound velocity] and different geometries 
(empty squares and dashed lines for the spatially homo- 
geneous case with periodic boundary conditions, filled 
circles and solid lines for the circular potential well). 

The e gaug0 -dependence in that figure allows to estimate 
the interval where the linear response approximation is 
reasonable, e.g. for e gaug0 / \JrakBTd < 0.08 the deviation 
due to non-linear effects is less than five percent. As ex- 
pected, the estimate (|35[) gives a more pessimistic bound 
around 0.04 for 7/(^5) = 0.4. From an experimental 
point of view, we expect that values of the gauge field 
amplitude as high as e gaU ge /(mk B T d ) 1/2 = 0.15 should 
be well achievable, see Table HI As shown in the figure, 
for e g augo/(wfcs Jd) 1 / 2 < 0.2, the effective normal frac- 
tion /„ empirically has, within the error bars, a linear 
dependence with e gauge - This could be expected from 
the fact that the first order correction to the linear re- 



sponse theory scales as e 



augc 



in the deposited energy 



AE, the term of order e gaugc in AE being zero due 



to the odd parity of the matter-wave coupling to the 
gauge field. Assuming this linear dependence of the effec- 
tive normal fraction, a fit was performed over the range 



(0.05) 2 < 



^/{mkBTd) < (0.2) , to extrapolate to 



fgaugc = 0, see the solid and dashed lines in FigSJ Note 
that this extrapolation procedure also eliminates the con- 
tribution to the effective normal fraction of the energy 
AEi/p deposited by the undesired light-shifts, see the end 
of Appendix [X] For larger values of the gauge field am- 
plitude e gauge , more serious non-linear effects set in in the 
effective normal fraction and it becomes less clear how to 



turn P y of the gas [4J,|45|]. The temperature dependence extrapolate them down to e gaug0 
of the normal fraction is shown in Fig|3] for increasing 
system sizes. The sudden variation around T/Td — 0.13 
becomes sharper and sharper as the system size is in- 
creased and should eventually correspond to a discontin- 
uous jump in the superfluid fraction at the BKT tran- 
sition [l|, 0, 5H . The slow decrease for larger values of 
T/Td is instead an artifact of the ultraviolet cut-off that 
has to be imposed to the classical field model in any di- 
mension d > 2. Indeed, the same decrease is visible also 



0. 



For the sake of completeness, it is important to note 
that for the weak gauge field amplitude that are required 
to be in the linear regime, or at least in the regime of 
the numerically suggested linear dependence of the ef- 
fective normal fraction with e gaugc , the deposited energy 
AE2 is at most a few percent of the total energy of the 
system, which may be experimentally challenging to mea- 
sure [43|. This value is however larger than the statisti- 
cal uncertainty of the energy in the canonical ensemble 
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FIG. 4: Classical field simulation of the deposited energy 
measurement scheme. Initial thermal equilibrium state at 
T/Td = 0.1 with g — 0.1 h 2 /m, and a system size L/£ ~ 63 
(the number of grid points is M = 64 2 , and n rca i = 1000 in- 
dependent realizations of the classical field are used). Gauge 
field modulation wavevector q ~ (l/£) e x . The real time evo- 
lution is followed during r = 3/7. The complete form Eq. (|21[) 
of the gauge field is considered, with ro located in the center 
of the system, and a beam waist w = 30£. For each value 
of the gauge field intensity e gaU ge of Eq. (|25[) , two calculations 
are performed to extract the energy change AE2 due to the 
spatially modulated component, from which an approximate 
value of the normal fraction /„ is obtained (see text): the first 
calculation with |f2+ | 2 (0 + )/|r2 c | 2 = e gaug e/(7ifc c ) and = 0, 
the second one with |n±| 2 (0+)/|n c | 2 = e gaugc /(2ftfc c ). Empty 
squares with error bars: spatially homogeneous system of size 
L (with periodic boundary conditions), which corresponds to 
N ~ 40000 atoms. Filled circles with error bars: system in 
a circular potential well (see text), with a total atom num- 
ber N ~ 30000 adjusted to have the same central density p, 
hence same degeneracy temperature Td and healing length £ 
as in the homogeneous case. The suggested experimental val- 
ues of Table U correspond to the circled point indicated by 
the oblique arrow and require the measurement of a relative 
energy change of ~ 0.7%. Thin lines: quadratic fit (with 
no linear term, see text) of the effective /„ as a function 
of £gauge ; over the interval 0.05 ^ fgauge /(mk B T d ) 1/2 < 0.2 
(dashed lines for the spatially homogeneous system, solid lines 
for the system in the potential well). Black, red and blue lines 
and points (from bottom to top) correspond to an excitation 
sequence with ^y/(c 3 q) = 0.4,0.2 and 0.1 respectively, with 
c s = (pg/m) 1 / 2 the Bogoliubov sound velocity. The extrap- 
olated zero-e gaug e values of the effective normal fraction are 
indicated on the figure. The region between the horizontal 
dotted lines indicates the confidence interval of the thermody- 
namic prediction shown in Fig(3] The residual 14% deviation 
of f n from the thermodynamic value (after linear extrapola- 
tion to 7 = 0) is a finite q£ effect (see text). 



with n rea i = 1000 realizations. As explained in the Ap- 
pendix [Ej such a large number of realizations turned out 
to be necessary in the small e gaU gc regime to compen- 
sate the fluctuations of the current j(t = 0) at the initial 
time. The resulting statistical error for the effective nor- 
mal fraction is about 2% for the parameters of Table [I] 



In the numerical results of Fig|4l this issue was circum- 
vented by a simulation trick so that n rea i = 1000 is a suf- 
ficiently large number of realizations to keep error bars 
below 2% in the limit e g auge — > [48|. It is also shown 
in the Appendix [E] that a measurement of the variance 
of the deposited energy (rather than its mean value) in 
principle also allows to access the normal fraction, if one 
knows experimentally the initial energy of the gas. 

After linear extrapolation to 7 = from the lowest 
two values of 7 in Fig [3] which gives /„ ~ 0.176, the 
residual disagreement of f n with the thermodynamic re- 
sult indicated by the dotted lines is about 14%, of the 
same order as the finite q correction predicted by Bo- 
goliubov theory, see Figj2^,. We have thus performed 
simulations for q£ = 1/2, for a spatially homogeneous 
system having the same parameters as in Fig|4] and for 
a gauge field switch-off rate j/(c s q) = 0.2 (not shown). 
It is found that the range of linearity of the effective 
normal fraction with £g aU gc is of the same order, up to 
egaugcAwA^Id) 1 / 2 ~ 0.27. Repeating the zero-e gaug e ex- 
trapolation procedure as in FiglU we find for q£ = 1/2 an 
effective normal fraction of 0.164(3), rather than 0.146(2) 
for q£ = 1. Assuming a quadratic dependence of the ef- 
fective normal fraction with q at fixed 7/(c s g), as pre- 
dicted by Bogoliubov theory, see Figl2jD, this leads to an 
additive correction to f n close to 0.024. We finally reach 
/„ ~ 0.20, which is within the statistical error bars of 
the thermodynamic value [5fl | . On a smaller system with 
M = 32 2 modes, we have reached a similar conclusion, 
taking even smaller values of g£ and 7 /{c s q) 1 see the blue 
circle in Fig[3] and the figure caption. The deviation ob- 
served in Fig[4] is therefore not a systematic error of the 
proposed method. 

As a final check, we have tested the locality of the 
proposed measurement scheme by performing the simu- 
lation for two different geometries. Empty squares and 
dashed lines in Fig|4] correspond to a spatially homoge- 
neous system with periodic boundary conditions, while 
the filled disks and solid lines correspond to a system 
trapped in a circular well with steep walls of the form 
U(r) = (k B T{tanh[(r - L/2)/(£/2)] + 1}, where the nu- 
merical parameter £ determines the energy cut-off in the 
classical field model 39] . The probed region is at the cen- 
ter of the potential well, ro = 0. Even in the non-linear 
regime, where linear response theory fails, the effective 
normal fractions are almost the same in both geometries. 



V. OPTICAL MEASUREMENT 

The proposal that we have illustrated in the previous 
section was based on the measurement of atomic quanti- 
ties, namely the deposited energy in the atomic cloud at 
the end of the gauge field sequence. The present section 
is devoted to the presentation and the characterization 
of an alternative, all-optical route to measure the normal 
fraction /„: information on the response of the atomic 
cloud to the gauge field can be retrieved from the trans- 
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2D density 


pA 2 = 9 


p = 14pm 2 


degeneracy 
temperature 


k B T d = 2irti 2 p/m 


T d = 500 nK 


temperature 


T = Q.lT d 


T = 50 nK 


2D interaction 
constant 


g = mg/h 2 


5 = 0.1 


transverse 
confinement 


hw z = 0.23 h 2 k 2 c /m 
= 0.16fc B T rf 


ft = 1.65 kHz 


healing length 


i = (P9)- 1/2 


£ = 0.84/im 




reduced gauge 
field amplitude 


£gaUgC - (mfc B T d )l/2 


tgaugc — 0.15 


probe beam 
Rabi frequencies 


2\n c \' J 


0.09 


gauge field 
switch-off rate 


7 = 0.2 c s q 
with q = \/( t 


1/7 = 4.8 ms 
q = 0.15/Cc 



level schemes of Fig|6] 


first choice 


second choice 


level |o) 


|F = 1,-1) 


|F = 2,-2) 


level \b) 


\F = 2, -2) 


1^ = 1,-1) 


level |e) 


|F' = 2,-1) 


|f = 2,-2) 


squared coupling Rabi 
frequency |f2 c | 2 /r 2 


0.21 


0.5 


minimum detuning 5/T 


1 


1.5 


fluorescence probability 
per atom Pfl uo 


0.22 


0.045 


spurious deposited energy 
AEu/AE 2 (for /„ = 0.2) 


33 


0.16 



TABLE I: Suggested values of the physical parameters for an 
experimental measurement of the normal fraction of a two- 
dimensional Bose gas of 87 Rb using an artificial gauge field 
produced by laser (coupling and probe beam) excitation on 
the Dl line with an optical wavelength A c = 795 nm. The first 
block characterizes the thermal equilibrium of the gas. The 
second block determines the gauge field. The third block deals 
with the issues of spontaneous emission and spurious light- 
shift, for two atomic level schemes: the "first choice" is Fig[6ji, 
the "second choice" is FigdJp. For this second choice, the 
indicated value |f2 c | 2 /r 2 = 0.5 corresponds to a compromise 
between minimization of the fluorescence and of the spurious 
light-shift, the truly minimal Pfl uo being < 0.01, sec (IA18I I. 
Note that Pfl uo and AEu / AE2 are basically unchanged if one 
takes q£ = 1/2 and j/(c s q) — 0.4. The three-dimensional 
scattering length 0,30 — 100 Bohr radii is related to the two- 
dimensional coupling constant by Eg . (|2Ufl . A useful relation 
is mT/(hk 2 ) ~ 792. c a = (pg/m) 1 ^ 2 is the Bogoliubov sound 
velocity. 



mitted probe beams once they have crossed the atomic 
cloud. Recent works have in fact pointed out that the 
strong frequency-dependence of the dielectric constant of 
an optically dressed medium in the electromagnctically 
induced transparency (EIT), already used experimentally 
to strongly reduce the light group velocity (5lTj53| . can 



be exploited for vclocimctry experiments: Information 
on the current profile of an atomic cloud was predicted 
to be imprinted onto the phase of the transmitted probe 
beam [Hj, HH • 

In the present case, the gas is illuminated by probe 
light of angular frequency oj p and coupling light of an- 
gular frequency w c . As in Scc lIVI thc coupling light is a 
plane wave propagating along y axis. In the subsection 
IV Al where it is proposed to measure the transmitted 
mean electric field amplitude, the probe light consists, 
as in Sec lIVl of two beams with wavevectors at a 
small angle with the z axis, as sketched in Fig{TJ dif- 
ferently from [Hil, [|55[ , the matter- wave current pattern is 
generated by the same laser beams that are then used for 
probing. The configuration is slightly different in the sub- 
section [VB] where it is proposed to measure the pattern 
of the scattered light intensity: in this case, the probe 
light consists of a single beam with a wavevector k p di- 
rected along the z axis, as sketched in FigJSJ Even though 
no spatially modulated matter-wave current is generated 
by the laser beams, still one can extract information on 
/„ from the zero-mean matter- wave currents induced by 
thermal fluctuations in the gas. 

A. Extracting f n from the amplitude of 
transmitted light 

The transmission and reflection of probe light from the 
two-dimensional atomic cloud can be described in terms 
of Maxwell equations. In particular, the dipole polariza- 
tion of the atoms provides a source term for the probe 
electric field E p at angular frequency uj p : for the positive 
frequency parts, one has in the paraxial approximation 
with respect to the z axis: 

(A + k 2 p )8 p = -^V P , (40) 

where k p = u p /c and the Laplacian operator is three- 
dimensional. In the considered A atomic configuration, 
see FigHJ one may think that the atoms occupy in each 
point of space the non-coupled state \NC(r)}, leading to 
a vanishing mean atomic polarization V p . This is ac- 
tually not exactly the case: Due to the atomic motion, 
the atomic internal state does not follow adiabatically 
and there exists a small coupling between the position 
dependent coupled state |C(r)) and non-coupled state 
\NC(r)}, the so-called motional coupling [30l ]. This leads 
to a small non adiabatic atomic polarization, that we now 
evaluate. 

Within a perturbative picture, we simply need to calcu- 
late the mean atomic polarization V p induced by the un- 
perturbed laser fields. In terms of the three-dimensional 
atomic field operators, this reads 

V p (r) = d ae (& a (r)* e (r)). (41) 

With the usual adiabatic elimination of state |e), as 
shown in the Appendix [B] the atomic field in the ex- 
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cited state can be written in terms of the atomic field 
operator \3D in the coupled |C) state, 

X3D (r, *) = (C(r, t)\a) *„(r, t) + (C(r, t)\b) * 6 (r, t) (42) 



as 



0, 



2(<5 + zr/2) 



[l + |fl P MP] 1/2 X3D+r 1/2 B, (43) 



From the explicit form of the noise term B e given in the 
Appendix [B] it is immediate to sec that it gives a zero 
contribution to the mean in Eq. (|41[) . Since the atoms 
arc mostly in the uncoupled \NC) state, we can approx- 
imate the atomic field \l/ a in the |a) state by its \NC) 
component. To lowest order in £l p /Q c we then have 



V P {v) 



2{5 + iT/2) 



(^L«X3D(r)>. (44) 



The next step is to perturbativcly evaluate the field X3D 
that is created in the coupled |C) state by the motional 
coupling [3(J of |C) to \NC). The details of the procedure 
are given in the Appendix iBl To first order in Cl p /Cl c 
and for the magic choice Eq. ([23| of A c , one obtains after 
adiabatic elimination 



y Vto ' v rf + r / ^ (45) 



As the proposal in [55| addressed a pre-existing current 
profile and the weak probe beam induced a vanishingly 
small gauge field, the difference was irrelevant in that 
case. Here, on the contrary, the mean current is itself 
proportional to the gauge field so that the difference be- 
tween the two operators really matters. Second, the ex- 
pression (|46j) contains an extra term proportional to the 
average density gradient. In contrast to the first term, 
this one is purely imaginary. As a result, it only affects 
the intensity of the transmitted light via a combination 
of absorption and/or amplification effects. In particular, 
as it does not induce any phase shift on the light, it does 
not interfere with the proposal of (55l |. 

In order to calculate the modification 5£ p induced by 
the atoms on the transmitted electric field of the probe, 
one has to insert the polarization (|46|) as a source term 
into the Maxwell equation (|4T)|) . Within a standard ap- 
proximation, we can neglect diffraction effects stemming 
from the in-plane part of the Laplace operator in (|40[) 
and integrate the z dependence across the atomic cloud. 
Taking into account the appropriate boundary conditions 
for S£ p , this leads to the expression 



6£ p ( X ,y,z)~e ik **^ 
zeo 



dz'e- lk " z VJx,y,z') (48) 



An explicit expression for the noise term B x is given in 
the Appendix [B] again, the noise term has a zero expec- 
tation value and does not contribute to the optical po- 
larization. The final form of the optical polarization in 
terms of the three-dimensional atomic density and cur- 
rent operators (TTT|) and (fTS")) reads, to zeroth order in 
q/k p : 



Vpir) _ x 



h\Sl r \ 2 



(j 3 D(r)) 



2im 



V(ji 3 £>(r)> 



£ p (r), (46) 



where the Rabi frequency of the probe beam has been 
eliminated in favor of the electric field using the defi- 
nition —d ea £ p = Ml p /2 and the detuning 5 has disap- 
peared from the formula. The first term proportional to 
the atomic current operator has a simple semi-classical 
interpretation in terms of the reduced group velocity in 
the Electromagnetically Induced Transparency regime, 
as anticipated in [5~I1. [55j: in this regime, the refractive 
index strongly depends on the Raman detuning, which in 
turn depends on the atomic speed because of the Dopplcr 
effect. 

However, the expression (|46p differs from the semi- 
classical one that was used in |55| in two ways. First, 
the current operator in Eq. (|4"6"f differs from the physical 
current of atoms by the gauge field 



jphys — J3.D n 3D A 3D . 

m 



(47) 



for the transmitted field in the z > region. In order 
for the approximation to be accurate, z has to be much 
larger than the thickness a^ of the atomic pancake, but 
at the same time much smaller than the diffraction length 
k p /q 2 , where q is the characteristic wavevector of the in- 
plane modulation of the atomic density and current. 

Along z, the atomic field varies as the harmonic oscil- 
lator ground state wavefunction 4>o(z), sec Eq. (fT9"|) . Per- 
forming the integral over z', this gives the final expression 
for the variation of the transmitted field 



S£ P (r) 



he \ fl c 



(jW) + ^-V(n(r)) 



£ P (r) 



(49) 

in terms of the two-dimensional density n(r) and current 
j(r) operators. The first contribution proportional to the 
atomic current gives a phase shift, while the second con- 
tribution proportional to the atomic density gradient is 
responsible for absorption and amplification of the probe 
beam. 

The atomic current profile created by the gauge field 
is evaluated using the linear response formulas (|3]) and 
((4]) as discussed in detail in the previous sections. The 
gauge field is assumed to be switched on slowly enough 
as compared to the characteristic frequencies of all the 
excitation modes of the gas at wavevector q. Within the 
linear response regime, the contribution to the current 
due to the spatially modulated gauge field at q may be 
isolated by a suitable combination of measurements with 
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different values of ft p , which gives 

0)2 (r) = (qcosa + f n e z x qsina) 

mq\\l c \ z 

x (n+n p * e jq - r + ci+*q~ e^ r ) e -*\*-*°\ 2 1™\ ( 50 ) 

where r is now in the xy plane and a is the oriented 
angle that k c makes with q [3y] . Inserting this expression 
into (|49|) and recalling the form of the incident probe 
field, one can extract the phase shift experienced by the 
central part of the probe beams after crossing the atomic 
pancake, 



= BY [cos 2 a + f n sin 2 a], (51) 

where we have assumed k p ~ k c . The total decay rate 
of the e state by spontaneous emission is indicated by 
r and B is the branching ratio for the decay to the a 
state, so that BY = \d ae \ 2 k^/ (Sirheo). As we have already 
mentioned above, the density gradient term in (|49[) only 
introduces an intensity modulation and is not responsible 
for any phase shift. 

From a nonlinear optics point of view, the phase shift 
(|51|) can be interpreted as arising from a X*- 3 -* optical non- 
linearity of opto-mcchanical origin similar to the one that 
was demonstrated in the experiment [56|: the nonlinear 
modulation of the optical response of the atoms is deter- 
mined by the mechanical distortion of the cloud by the 
optical forces. 

Inserting into (fS"Tj) the values of Table fl] for the 87 Rb 
case, for the first choice of level scheme reported in the 
Appendix [A] one has a branching ratio B = 1/4 and one 
finds a small, yet appreciable phase shift on the order of 



B. Current fluctuations and the angular 
distribution of scattered light 

All the calculations presented in the previous subsec- 
tion aimed at evaluating the expectation value of the 
transmitted field amplitude. At this level of the theory 
we were allowed to describe the probe beam as a coher- 
ent, classical field and we could neglect the fluctuations 
around the expectation value of both the light field ampli- 
tude and the atomic current and density operators. The 
formalism can be straightforwardly extended to quantum 
optical fields so to include the fluctuations of the atomic 
density and current. This is crucial when one aims at 
investigating the spontaneous scattering of light off the 
current fluctuations in the atomic gas. In this subsec- 
tion, we shall in particular show how information on the 
normal fraction of the gas can be inferred from the an- 
gular distribution of scattered light. We shall make the 
approximation of replacing temporal derivatives of the 
electromagnetic field dt£ by ~ick p £ in Maxwell's equa- 
tion. In particular, this misses retardation effects in the 
expression of the scattered fields in terms of the atomic 
dipoles, which is accurate since the system size is much 
smaller than c/Y. 

jz 

--■ r- p 



x 



A0: 



choice 1 



6-10 4 x [cos 2 a + f n sin 2 a]. 



(52) 



For the second choice of level scheme in the Appendix 
IA1 the branching ratio is slightly larger, B = 1/3, but 
for the compromise choice Eq. (|A25[ ) |il c | 2 /r 2 is larger so 
that one finds a smaller phase shift 



A<^ hoico 2 



3 ■ 10~ 4 x [cos 2 a + f n sin 2 a]. 



(53) 



In addition to the phase shift of the transmitted beam 
that we have discussed so far, Bragg diffraction on the 
spatially modulated current profile produces a pair of ad- 
ditional beams of in-plane wavevector respectively ±3q/2 
via a sort of four- wave mixing process. The relative in- 
tensity of these beams as compared to the incident probe 
beams is on the order of |A(/)2| 2 . For transverse gauge 
fields such that q • k c = the contribution of the induced 
density gradient term of (|49[) vanishes by symmetry. In 
the case of longitudinal gauge fields, the relative correc- 
tion is on the order of q£. 



FIG. 5: Scheme of the scattering geometry under exami- 
nation. Probe light is incident at wavevector k p = k p e z 
and the scattered light is collected at a wavevector k sc = 



(k p 



2^1/2, 



+ Q 



We consider the geometry sketched in Fig(5j a single 
Gaussian probe beam is incident onto the atoms with a 
wavevector k p exactly orthogonal to the atomic plane, 
waist w centered at ro and a weak peak amplitude £ p . 
Inserting this form into (|49[) and taking the Fourier trans- 
form along the xy plane, one obtains the following opera- 
tor equation for the scattered field component at in-plane 
wavevector Q 57], 



£(Q,z) 



2z kp \da 
e h |O c 



ttw 2 £° 



d 2 q 

(27T) 2 



, fi(Q-q) 
jQ - q + 2m ^ 



Jk z (Q) z 



(54) 



Here, jq and nq are the spatial Fourier transforms of 
the two-dimensional current j(r) and density n(r) opera- 



13 



tors; the Fourier transform of a product of two functions 
has been rewritten in terms of the convolution of their 
Fourier transforms. The z component of the propaga- 
tion wavevector is determined by the photon dispersion 
as k z (Q) = (fc 2 - Q 2 ) 1/2 - 

The intensity of the scattered light at in plane wavevec- 
tor Q is quantified by 



(£t(Q)£(Q)) = 
d 2 q f d 2 q' 



2nk p \d ae \ 2 W 2 \£° 



(2tt) 2 J (2tt) 2 
x Ik 



jQ-q 



HQ - q) 

2m 



, ft(Q-q') 

jQq ' + 2m nQ " q '. 

(q-q')-rn -(g 2 +g' 2 



x e 



V4 



(55) 



Since the system size is much larger than the waist w of 
the probe beam, we can for simplicity assume an effective 
translational symmetry along the xy plane. As a con- 
sequence, the correlation function that appears in (|55j) 
has a delta-distribution shape around equal wavevectors 



Q q = Q — q' (see e.g. the next Eq.([56|l). 

In contrast to the schemes proposed in the previous 
sections, where the duration 1 of the experiment had to 
be at least on the order of l/(c s g), where c s = {pg/m) 1 / 2 
is the Bogoliubov sound velocity, the light scattering ex- 
periment discussed here can be performed on a much 
faster time scale, only limited by the characteristic rate 
r' of the internal atomic evolution time, Eq. (|B7[) . As 
a result, the experiment can be performed in the small 
wavevector region Q£ <C 1 where the contribution to |55|) 
of the terms involving the density fluctuations uq is neg- 
ligible [58}. Of course, efficient isolation of the scattered 
light from the incident beam requires that the scatter- 
ing angle 9 SC ~ Q/k p be much larger than the diffraction 
cone of the probe beam, i.e. Q ^> l/w. 

The instantaneous correlation function of the current 
in the y direction parallel to k c can be evaluated applying 
the fluctuation-dissipation relation (|E5[) to the current 
operator j q in an infinite space geometry. This gives 

(jQ, y JQ',v) = (27r) 2 S 2 (Q-Q')k B TR C [ X yy(Q',Lo = 0)] 



Q-s-0 



(27r) 2 ^(Q- Q 



,.k B Tp(r ) 



x [cos 2 <f) sc + f n sin 2 SC ], (56) 



where <j) sc is now the azimuthal angle between k c and 
Q. Inserting this expression into (|55[) and taking the 
thermodynamic limit, one gets to the final expression for 
the scattered intensity in the momentum Q-space (57| . 



e n\n c \ 2 

k B Tp(r ) 



m 



2tiV \£ 



(I I 2 



cos 2 4> sc + f n sin 2 



(57) 



To estimate the relative intensity of scattered light, it 
is useful to rewrite the expression (|57[) for the momen- 
tum space intensity (£q£q) in terms of physically more 
transparent quantities such as the angular distribution 
I(0 sc ,4>sc) of scattered intensity. For small scattering 
angles |0 SC | <C 1, the infinitesimal solid angle and mo- 
mentum space volume elements are related by dQ = 
sin 9 SC dd sc d(j) sc ~ 9 SC dd sc d<j) sc ~ d 2 Q/k 2 , so that 



(27T) 



(58) 



This immediately leads to the final expression for the 
angular distribution of the scattering intensity [HtJ 



^sc (^sc i 0sc ) 



kp k c \d ae 



k B Tp(r ) 



ne h \Q C \ 2 

[cos 2 SC + f n sin 2 SC ] (59) 
in units of the incident intensity, 



Itac = Jd 2 r\£ p (r)\ 2 ^^f\£° p \ 2 . 



(60) 



From this expression, it is immediate to see that informa- 
tion on the normal fraction of the gas can be retrieved 
from the azimuthal dependence of the scattered inten- 
sity. In terms of the total spontaneous emission decay 
rate T of the e state and the |e) — > \a) branching ratio B 
(we recall that BT = |d ae | 2 ^/(37rSe )), Eq.([59"]) can be 
rewritten in the more transparent form [57| . 



9k 2 B TT d B 2 T 2 

2wh 2 in. 



-'-sc \^sc i Y sc J 



[cos 2 0sc + /„ sin 2 



(61) 



To estimate the feasibility of the proposed light scat- 
tering experiment, we now derive an upper bound on 
the number of useful scattered photons in a single shot 
of duration r. Calculating the Poynting's vector of the 
probe beam, and using (|60|). we find an incident flux of 
probe photons $ inc = (k p w) 2 \fl^\ 2 / (12BT) . Integrating 
in (|59|) the term proportional to /„ over solid angles in 
the cone 9 SC < l/(k p £), we obtain the flux 4>^ e of use- 
ful scattered photons. As a maximal duration, we take 
t = l/r{5°"~ ad where the fluorescence rate of the atoms 
due to motional coupling between the non-coupled and 
the coupled states is given by Eq. (lA2l) (with 2|£7+| 2 re- 
placed here with |fip| 2 ). The number of single shot useful 
scattered photons is thus bounded by 



N, 



ph 



16 hjj z 



(keif 



(62) 



Remarkably the Rabi frequencies fl c and f2° have can- 
celed out in the ratio of the scattered flux to the fluores- 
cence rate. One recognizes in the right-hand side of (f6"2")l 
the effective mean number of atoms 7V° t ff = irpw 2 /„/4 in 
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the normal component illuminated by the probe beam, 
as in Eq. ([2"B)l . There is however a severe geometrical re- 
duction factor, l/(/c c £) 2 , due to the small aperture of the 
useful scattering cone. For the parameters of Table U 
with B = 1/3, and taking a waist w = 30£ and /„ = 0.2 
as in FigH we find iVff ~ 1400, l/(fc c £) ~ 0.15, which 
leads to N^ e < 5. This remains accessible to current 
quantum optics experiments. For ksT/huz fixed, the 
upper bound in Eq. (|62|) scales as p 2 , since l/£ 2 scales as 
p, so that larger values of photon numbers for a given 
waist may be obtained by increasing the density p of the 
two-dimensional Bosc gas. 

VI. CONCLUSIONS AND PERSPECTIVES 

In this paper we have proposed and validated two 
methods to measure the superfluid fraction of a quantum 
fluid of ultracold atoms. The idea is to apply an artificial 
gauge field to the atoms with spatial oscillations within 
a localized envelope and to detect, within the linear re- 
sponse regime, the matter-wave current pattern that is 
generated in the fluid. This can be done in a mechan- 
ical way by measuring the energy that is deposited in 
the fluid at the end of a gauge field sequence. This can 
also be done in an all-optical way by observing the phase 
shift experienced by the same laser beams that are used 
to generate the artificial gauge field or the angular pat- 
tern of scattered light. We have shown that, by a careful 
choice of the parameters and of the atomic level scheme 
for 87 Rb, two experimental obstacles, the spontaneous 
emission and (in the case of the mechanical method) un- 
desired light-shifts, can be put to an acceptable level. 
The experimental challenge remains in the required high 
sensitivity of the measurements, that is a detection of 
small energy changes in the mechanical method, and a 
detection of small optical phase shifts or small photon 
numbers in the all-optical methods. 

The interest of the proposed methods is twofold: they 
do not require that the gas reaches thermal equilibrium 
in presence of the gauge field, and furthermore they give 
the possibility of reconstructing in a local way the spa- 
tial profile of the superfluid fraction of a trapped gas, 
independently from the presence or the absence of a Bosc- 
Einstcin condensate. This last feature is attractive in the 
study of the Bcrezinskii-Kostcrlitz-Thouless transition to 
a superfluid state in two-dimensional Bose gases and of 
the superfluidity properties of Bose gases in disordered 
environments. It would also be interesting to extend the 
method to the study of superfluidity in multi-component 
atomic fcrmionic gases, which may require identification 
of suitable level schemes. 
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Appendix A: Experimental issues 

In this Appendix we review some issues that may hin- 
der an experimental implementation of our proposal. Our 
attention will be concentrated onto the most relevant case 
of 87 Rb atoms considered in the experiment of [HiJ |60| 
and, in a two-dimensional context, in The level struc- 
ture of this atomic species is sketched in Fig|5J A possibly 
important advantage of this atom in view of an experi- 
mental implementation of the present proposal is that 
the singlet and triplet scattering lengths of ground state 
atoms are equal within a few percent |6lj . We thus ex- 
pect that all scattering lengths between arbitrary F = 1 
or F = 2 sublevels of the ground state have all almost the 
same values (62|, which leads to g aa ~ g a b and therefore 
to a suppression of the spatio-temporal variation of the 
effective interaction constant g%D (r,i) defined in (fl~8)). In 
the deposited energy method proposed in Scc lIVl this 
is important to reduce the emission of phonons in the 
atomic gas by the temporal modulation of the interac- 
tion constant. In the optical detection scheme of SccfVl 
this is also important to suppress the contribution of the 
interaction term to the field x^d and then to the opti- 
cal polarization. Other atomic species such as Yb [18[ 
or metastable He I63l . l64i and/or different laser beam 
configurations [l9|, [2(J |26| are expected to be useful for 
other purposes, e.g. to suppress spontaneous emission 
and/or generate artificial gauge fields with different ge- 
ometries [65| . 

Two possible choices for the three states |a),|6),|e) 
forming the A system on the Dl line of 87 Rb are con- 
sidered, as sketched in the two panels of FiglfJ] For each 
choice, we determine the undesired effects (spontaneous 
emission, light-shifts, Raman leaks) stemming from de- 
viations from the perfect adiabatic following of the non- 
coupled state by the moving atoms and from optical tran- 
sitions to other levels not included in the A system. An 
eye will also be kept on trying to maximize the |e) to 
| a) branching ratio so as to reinforce the optical signal of 
Sec|V] To minimize spontaneous emission within the low 
saturation regime, we shall allow for a small detuning 6 
of both the probe and coupling beam carrier frequencies 
from the \a) — > \e) and \b) — > |e) transitions, respectively. 
The Raman detuning of the two beams is taken in a way 
to always fulfill the magic Raman condition (I2U1) . 

a. First choice 

The coupling beam propagates along the y axis and 
is taken with a cr+ polarization with respect to the y 
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FIG. 6: Scheme of the internal levels of 87 Rb atom involved in the Di transition (J = 1/2 — >■ J' = 1/2). The hyperfine splitting 
of the ground (excited) state is Ahf, g = 2-7T x 6.834 GHz (Ahf, c — 2n x 814 MHz). The natural linewidth of the excited state 
is r = 27T x 5.75 MHz. Thin (thick) lines indicate transitions that are induced by the probe (coupling) beams. Solid (dashed) 
lines indicate the desired (main undesired) transitions. The relevant dipole matrix elements are shown in units of the reduced 
dipole element (J = l/2||er|| J' = 1/2) of the Dl line, in boldface for the desired transitions and in italic for the undesired 
ones. Our first choice for the A system \a), e), \b) states is shown in the left (a) panel, with \a) = \F = I.tuf = —1), |e) = 
\F' = 2,m F i = —1), 1 6) = \F = 2,m F = —2). The second choice is shown in the right (b) panel, with \a) = \F = 2,iuf = —2), 
\e) = \F' = 2, m F i = —2), |6) = \F = 1, tuf = —1). Note that we have taken the y axis as the quantization axis of angular 
momenta. 



quantization axis. The probe beam is taken as linearly 
polarized along y. For the three atomic levels forming 
the A system, we take \a) = \F = l,mp = —1), \b) = 
\F = 2,m F = -2) and |e) = \F' = 2,m F > = -1) [66|. 
This scheme of levels and lasers is illustrated in the left 
panel (a) of FiglHl To estimate the importance of the 
non-adiabatic coupling between the \NC) and the \C) 
states due to the atomic motion, we can evaluate the ra- 
tio pc/pnc of the two-dimensional densities in the two 
states close to the center of the laser spot. This is done 
using the explicit formula (|45|) for the field in the \C) 
state [69|. Using the fact that for quasi-2D samples the 
gradient is mostly along the harmonically trapped z di- 
rection and assuming n+ = Q~ , we obtain for k p ~ k c : 
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4(S + iT/2)k c n+ 
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TO 



Accuracy of the adiabatic approximation requires that 
this ratio is much smaller than unity. 

The finite population that is present in the |C) state as 
a consequence of non-perfect adiabaticity is responsible 
for the spontaneous emission of photons at a single atom 
rate: 



iiion-ad 
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where the fluorescence rate T' of the coupled state |C) is 
defined by Eq. (|TT|) and a spatial average has been per- 
formed. Remarkably, r^°"~ ad does not depend on the 
detuning <5. 



Other contributions to the fluorescence rate come from 
non-resonant excitation processes. Dominating among 
these are the excitation of the \a) state to the states \F' = 
1 or 2, to pi == 0) by the coupling beam at a rate 



^fluo — r 



in, 



(A3) 



and the excitation of the non-coupled state \NC) to 
\e') = \F' = l,mpi = —1) by the total probe plus cou- 
pling field with an effective Rabi frequency — 2n p /v / 3, 
which results in the fluorescence rate on the parasitic A' 
configuration \a) — > \e') — > \b): 
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In these expressions, we have taken into account the tab- 
ulated hyperfine dipole matrix elements of the various 
optical transitions, and we have introduced the hyperfine 
splittings Ahf, g and Ahf, e given for 87 Rb in the caption 
of FigUU Limiting ourselves to the most relevant regime 
where |n+/n c | 2 > 1/100, we see that in the present case 



of 87 Rb atoms, » rg uo . 

The total fluorescence rate can then be approximated 
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by a careful choice of the coupling beam Rabi frequency 



\n° pt \ 



A 
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(A6) 



Here, we have expressed lu z in terms of the three- 
dimensional scattering length a-^jj and the reduced two- 
dimensional coupling constant as given by Table U and 
Eq. (|20)l . Inserting the actual parameters of the 87 Rb 
atom, and taking g = 0.1, wc obtain |n° pt | 2 /r 2 ~ 0.21. 
It remains to adjust the detuning 5 to be in the weak 
saturation regime, 
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Let us take the same values of the gauge field se- 
quence as in FigHJ j/(c s q) = 0.2 and q£ = 1 and 



0.15{mk B T d ) 



1/2 



{pg/m) 1 / 2 being the Bo- 



goliubov sound velocity. This choice of e ga uge leads to 
\n+(0+)/n c \ 2 ~ 0.03(pA 2 ) 1/2 : as this quantity has to be 
much smaller than 1 in order for the gauge field descrip- 
tion of Sec lIIII to be valid, it is safe to impose pX 2 < 10. 
Integrating over the exponential switch-off ramp of Q p 
and eliminating g in terms of oj z and 0315, this gives for 
the total fluorescence probability per atom, 
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(A8) 

For 87 Rb with the choice pX 2 = 9, we obtain the not very 
impressive result, 



Pa* 



0.22. 



(A9) 



For the sake of completeness, it is important to note that 
for this choice, k B T — 0.1k B Td remains smaller than Huj z , 
so that the Bose gas retains a two-dimensional character. 

The existence of other atomic levels in addition to the 
ones strictly needed to create the gauge field is respon- 
sible not only for dissipative effects such as fluorescence, 
but also creates reactive effects such as a spatially and 
temporally-dependent light-shift of the non-coupled state 
\NC). Among the dominating processes, the parasitic A' 
scheme creates a modulated light-shift potential 
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ft|Q P (r)| : 
3A h f, c 
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A shift of the same order of magnitude arises from the 
coupling of the \a) state to the \F' = 2 or 1, mpi =0) by 
the coupling beam. 

An estimate of the energy deposited in the system by 
the U A term as compared to the one AE 2 due to the 
gauge field can be obtained with Bogoliubov theory: Us- 
ing Eq.dmll) with U Q = R(n+n-*)(0+)/(3A hf) e), v = 7 
and Q = q, one gets for Try <C e q <C pg: 



AE. 



JJA> 



la 



AE 2 



9/„(fc c c s ) 2 A 2 fi 



(All) 



For actual parameters, the energy change due to AU A 
turns out to be non-negligible. For the optimal value of 
the Rabi frequency f2° pt from Eq. (|A6j) . the ratio is 
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2hco z 
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For 

AE; 



87 Rb, one finds the discouraging result 
* / AE 2 — 600<?/( / oA 2 /„), which remains much 



larger than unity even for pX 2 = 9. 

Even if the deposited energy by the spurious potential 
U A is much larger than the desired one of the gauge field, 
a suitable extrapolation procedure may take advantage of 
the different dependence on the laser intensities to isolate 
the effect of the gauge potential. An alternative possibil- 
ity is to exploit the fact that the AE UA > contribution docs 
not depend on the direction of q: Within the regime of 
linear response, this contribution can therefore be elimi- 
nated by taking the difference of the energy changes for 
respectively longitudinal and transverse gauge fields. 

Another possible nuisance is the existence of stimu- 
lated Raman processes that may out couple the non- 
coupled \NC) state to atomic ground state sublevels |c) 
other than la) and \b), via excited state sublevels other 
than |e) [7PJ. One may however check that for the pro- 
posed scheme these leaky Raman couplings are detuned 
from resonance by a frequency amount at least Ay g in 
absolute value and therefore harmless. 



b. Second choice 

Another possible choice for 87 Rb atoms is to take a 
cr_ polarization for the coupling beam propagating along 
the y axis. The probe beam is again linearly polarized 
along y. The atomic levels forming the A system are now 
I a) = \F = 2,m F = -2), \b) = \F = l,m F = -1) and 
|e) = \F' = 2,m,F' = —2). The strong two-body losses 
that are generally experienced by the upper hyperfinc 
manifold of the ground state are here suppressed by the 
choice of a maximal mp state for \a): as collisions be- 
tween ultra-cold atoms mostly occur in the s-wave scat- 
tering channel, conservation of the sum of the mp's then 
prevents transition to the lower hyperfine manifold. 

The fluorescence rate per atom due to the motional 
coupling between \NC) and |C) is still given by r£°^ ad 
as defined in Eq. (|A2j) . As there is no longer any parasitic 
A' system, the fluorescence due to laser excitation of |o) 
or I b) to excited state sublevels other than |e) is now dom- 
inated by the transitions \b) —> \F' = 1 or 2,mp/ = — 1) 
due to the probe beam. Thanks to the reduced occu- 
pation probability ~ |17 p | 2 /|O c | 2 of sublevel \b) in the 
atomic state \NC) and to the larger hyperfine splitting 
Ahf )g of the ground state, the fluorescence rate is strongly 
suppressed. After spatial averaging it amounts to 
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Other fluorescence processes on the D2 line (e.g. the 
transition \a) — » |J' = 3/2, F' = 3,m F i = —3) excited 
by the coupling beam) are several orders of magnitude 
weaker than T^ no thanks to the huge fine structure split- 
ting of 2tt x 7 THz. 

Since the two terms in the sum Tg°"~ ad + Fg uo experi- 
ence different switch-off functions e -7 * and e~ 27 *, we in- 
tegrate over time to calculate the total fluorescence prob- 
ability: 



4 r|n+(o+)| 2 ^ 3r|n+(o+)| 4 
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After optimization over the coupling beam intensity, the 
minimal fluorescence probability 
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is obtained for a coupling beam Rabi frequency such that 



|n° pt | 4 
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Introducing the reduced quantities 7 = frj/(pg) and 
fgaugo = EgaugeAmfce'Td) 1 / 2 , and eliminating lo z in terms 



of g and a^n, we finally obtain 
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r fluo 
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For the parameters of Table U in particular p\ 2 = 9, one 
finds |fi° pt | 2 /r 2 ~ 5.5 so that a detuning \5\ > hT is 
required to remain in the weak saturation regime. This 
resulting probability of spontaneous emission per atom 
in the deposited energy measurement is very small, 



0.008. 
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As compared to the first one, this second choice then 
provides a strong reduction of the spontaneous emission 
rate by a factor almost 30. 

Another advantage of this second choice is that light- 
shift effects are potentially smaller thanks to the absence 
of the parasitic A' scheme. The probe beam on the \b) — s- 
\F' = 1 or 2,mpi = —1) transitions produces a light- 
shift which, after average in the \NC) state, leads to the 
spurious potential 
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The amount of energy that is deposited in the gas by 
this spurious potential can be estimated using twice 
Eq.(|Dl5|, first with U = -h\Sl+(t = 0+)| 4 /(|Q c | 2 A h f, g ), 
r) = 27, Q = q, and second with Uq four times smaller, 
j] = 2j, Q = 2q. Neglecting 2tvy with respect to e q , and 
taking = 1, we obtain the following estimate for the 
spurious deposited energy 
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For the coupling beam Rabi frequency (|A16[) minimizing 
spontaneous emission, the ratio of the energies deposited 
by the gauge field and the spurious potential amounts to: 
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For g = 0.1 and e ga uge = 0.15 and using the 87 Rb pa- 
rameters summarized in Table HI one finds the still quite 
unfortunate result 
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A possibility to overcome this difficulty and separate AE 2 
from AU P is to use the same strategy proposed to sep- 
arate AE 2 from AEi by exploiting the different depen- 
dence of the two quantities on the ratio |f2+/f2~|(0 + ). 
This can be combined with the choice of a compromise 
value of |r2 c | 2 /r 2 that allows to strongly suppress the 
light-shift potential without introducing a too large spon- 
taneous emission rate. 

To this purpose, we fix e ga ugc = 0.15, 7 = 0.2pg/fi, 
q£ = 1, g = 0.1 and we take as free parameters X = pA 2 
and Y = |f2 c | 2 /r 2 . Inserting the relevant parameters for 
87 Rb as in Table HI we obtain 

P«„o = 6 ' 86 :i°r 2 + 7.46 • 10- 4 F, (A23) 
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A reasonable compromise between the two competing ef- 
fects is to choose X = 9 and Y = 0.5, which corresponds 
to 



pA 2 = 9 and 



10, 



r 2 



= 0.5 



(A25) 



As a result, for the same parameters T/T^ = 0.1 and 
fn = 0.2 used in the classical field simulations of section 
IVl we obtain the quite encouraging values 



-Pfluo 
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We have checked that at the resulting temperature 
ksT/huj z = 0.62 the Bose gas retains a two-dimensional 
character and that the validity of the gauge field model 
of Sec lIIIl is guaranteed by the resulting probe beam Rabi 
frequency |fi+(0+)| 2 /|n c | 2 = 0.03A 1 / 2 = 0.09 < 1. 

To conclude, we remark that the energy deposited by 
the spurious light-shift is suppressed in the extrapolation 
procedure proposed in subsection IIVCI to eliminate ef- 
fects beyond the linear response regime to the gauge field, 
that is an extrapolation of the effective normal fraction 
that is linear in e 2 augo - Since AEup scales as |n+ (0 + )| 8 
with the probe beam Rabi frequency, it varies as £g auge 
with the amplitude of the gauge field, so it indeed con- 
tributes to the effective normal fraction as eg aug0 , which 



18 



extrapolates to zero. This works of course for e gaU ge small 
enough for the mechanical effect of the spurious light- 
shift to be treatable in the linear response regime. 



Appendix B: Dum-Olshanii theory for many-body 

systems 

In a seminal work [l7j , Dum and Olshanii have shown 
that an effective gauge field appears in the theoreti- 
cal description of a three-level atom interacting with a 
laser field on a A transition. Here we use the formalism 
of the Quantum Stochastic Differential Equations (see 
e.g. §8. 3. 2 of HH) to extend this idea to an interacting 
Bose gas in second quantized form. 

We start with the master equation for the density oper- 
ator a of the many-body system, assuming for simplicity 
that spontaneous emission corresponds to a net loss of 
atoms, 

±a=±[H,&] 
dt ih l J 



r J d 3 r * e (r)a#t(r) - l{$t(r)$ e (r),v} 



(Bl) 



where [, ] ({, }) stands for the commutator (anticommuta- 
tor). The Hamiltonian H is the sum of the single-particle 
kinetic and trapping terms, of the interaction terms (that 
we formally model as local Dirac-delta interactions), of 
the internal energy of the atomic excited state, and of 
the coupling terms of the atoms to the laser fields. 

Since the loss rate KT greatly exceeds the kinetic, trap- 
ping and interaction energies, we can neglect the external 
dynamics of the excited state and write 
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As previously defined, S is the common value of the 
detuning of the probe and coupling beams from the 
| a) — > \e) and \b) — > \e) transitions. 

In a Heisenberg picture for the open atomic system, 
the ground state atomic field operators ^ a =a.b satisfy 
the usual evolution equations ihd t ^> a = H]. On the 
other hand, conservation of the canonical commutation 
relations of the fields and of the Hcrmitian conjugation 
relation between if? a and \&t requires including a quan- 
tum Langevin term F e in the evolution equation for the 
excited state field ^p, 



Here, the quantum noise term F e is (S-correlated in 
position and time, e.g. [F e (r, t), F^(r', t')] = d(r — 
r')<5(t — t'), and we recall that the expectation value of 
normally-ordered products of noise operators vanish, e.g. 
(F% F e ) = 0, since the bath does not provide an incoming 
flux of e atoms. 

The only non-zero contributions to the commutator in 
Eq. (IB3|) originate from the excited state internal energy 
and from the atom-laser coupling term. This latter term 
can be expressed solely in terms of the atomic field op- 
erator x.3D in the coupled internal state |C), as defined 
in Eq.([4l?j). Along the lines of [7l|, we formally integrate 
[d t + (—iS + r/2)]*I' e = S neglecting a transient of dura- 
tion 1/r as 



r+ao 

*eM)= / dTe- { - u+T/2)T S{v,t-T). 
Jo 
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The Rabi frequencies f2 CiP and the atomic field X3D have 
a negligible variation during 1/r and may be replaced 
by their values at time t in the integrand. This leads to 
Eq. (|43"|) of the main text, where the noise term is defined 
as B e (r, t) = J +o ° dr e"(- j5 + r / 2 )-F e (r, t-r). 

As explained in section IIII| we are in a regime where 
the atoms are mostly in the non-coupled state and the 
field X3D in the coupled state is small and a perturba- 
tion expansion in powers of X3D can be performed. The 
gauge field formalism discussed in Scc lIIII for the evo- 
lution of the atomic field (f> 3 D in the non-coupled state 
[defined in Eq. ([T2"j) ] is already recovered at zeroth order 
in X3D- From this zeroth order approximation of 03_d, it 
is then easy to obtain the first order contribution to the 
field X3D that is required in Sec|V]to evaluate the optical 
polarization of the moving atoms. 

From Eq. (|12[) the equation of motion for c/j^o is 



(B5) 



d t 4>3D = ^[hn,H]+ ]T i> a at(NC\a). 



a—a,b 



By the very definition of non-coupled state, the excited 
state internal energy and the atom-laser coupling terms 
give an exactly vanishing contribution to the commuta- 
tor. In the other terms of the Hamiltonian as well as 
in the last sum in Eq. (|B5[) . we can perform the approx- 
imation ty a ~ (a\NC)<f>3D, which is accurate at zeroth 
order in X3D- After an integration by part and noting 
that (NC\d t \NC) and (NC\V\NC) are purely imaginary 
quantities, we find that up to this order fop follows a 
purely Hamiltonian evolution governed by Eq. (|13p . 
The equation of motion of X.3D has the form 



dtX3D = \[X3D,H] + Vc&iQa) 



(B6) 
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The commutator with the internal excited state energy 
term introduces a ^> e term, that we replace with Eq. (|43| : 
in this way, both a noise term and a complex, position 
dependent energy term —{i8' + V /2)x3D appear in the 
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equation. The real quantities 5' and V are given by 
Eq. (fTTj) and correspond to light-shift and damping ef- 
fects, respectively. 

Since HT' is much larger than the kinetic, trapping, in- 
teraction and recoil energies of the atoms, we can neglect 
these latter terms in the evolution equation of the coupled 
state, and only keep the coupling to </>3.d. This amounts 
to keeping in Eq. (|B6j) only the contributions to the ki- 
netic, trapping and interaction terms of the Hamiltonian 
H that contain one single factor x\d ano - an arbitrary 
number of ip^D and factors. In this way, we obtain 

d tX 3D * -(tS' + ^) X 3D+4>3D(C\[-d t + ^A]\NC) 

+ -V0 3D ■ (C\V\NC) + ^g$Ul D + T' 1 ' 2 F X (B7) 
m in, 

where we have introduced a complex position and time 
dependent coupling constant 

Q = (C\a)(a\NC)[(g aa - g ab )(NC\a) (a\NC) - a o b] 

(B8) 

with the convention gi> a = g ab . The noise term is de- 
fined by F x = -i\S + iT/2\(fl*/\n c \)B e . Its correlation 
properties are determined by the commutation relation 
[F x {r,t),F^{r',t')] ~ |w r r/a| ' e i{t-t')S-v\t-t'\/2 5 ^ _ r ,y 
Since we are working in a low saturation regime in which 
r ^> F', the time dependent factor in front of <5(r — r') 
may be replaced with a Dirac of t 
practice a spatio-temporally delta-correlated noise. 

The last step is to expand Eq. ([B7[) to first order in 
f2 p /f2 c . We also limit ourselves to zeroth order in q/k p 
and in l/(k p w), and we neglect the temporal derivative 
of the switch-off function f(t). Then (i) for g aa ~ g ab , 
the fourth contribution in the right-hand side of Eq. (|B7p 
vanishes [72|], and (ii) for the magic choice Eq. (f23|) , the 
second contribution in the right-hand side of Eq. (|B7j) 
vanishes. With the same adiabatic elimination technique 
adopted for ^> e and taking into account the fact that 5' 
and r' vary very slowly on the scale of 1/r', we are fi- 
nally led to the final equation Eq. (|45| with a noise term 
defined by B x (r,t) = J + °° drer^ 6 ' + v ' / 2 ^^ T F x {v, t - r). 

We complete the discussion by giving the back-action 
of the field \3D on the field 4>3D, a back-action that was 
already considered for a specific single atom geometry in 
[30| . The linear coupling of 3 £> to X3D originates from 
terms in the Hamiltonian that are linear in X3D, leading 
to 
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Expression of the back-action solely in terms of 4>3D and 
noise operators is obtained by replacing X3D in the result- 
ing equations of motion with its adiabatic approximation 



derived from Eq. ()B7p . This leads in general to a lengthy 
formula. For simplicity, we give the result for g aa = g a b, 
to leading order in f2 p /f2 c , we also neglect the contribu- 
tion to (C\dt\NC) of the time-derivative of the switch-off 
function f(t), and we use the specific form Q p /Q c con- 
sidered in this paper, see Eqs. (|6l7p . restricting to zeroth 
order in q/k p and l/(k p w) so that 
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After reduction to the xy plane, the deterministic term 
gives rise to two corrections to the evolution of 3 £>: (i) 
a complex position dependent energy shift, 

h{5" ,T"/2) = -2(6 + ff/2K^M, (Bll) 

m |SZ C | 4 

and (ii) a complex correction to the mass along y, 5m v = 
8hk 2 \n p \ 2 (S + iT/2)/\Q c \ 4 . The quantity hS" is the light- 
shift potential experienced by the non-coupled bidimen- 
sional field. The spatial average of the fluorescence rate 
T" of the non-coupled field coincides with the r£™- ad flu- 
orescence rate previously discussed in (|A2|) . as it should 
be. For the parameters of Table Q] the reactive correc- 
tions 6" and Sm y are small provided that the detuning is 
not too large, \5/T\ < 5. For instance, an estimate for the 
undesired energy deposited by the HS" potential can be 
obtained from the Bogoliubov theory for a homogeneous 
system, see Eq. (|D15p . leading to 
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where AE 2 is the desired deposited energy giving access 
to the normal fraction /„. 



Appendix C: Derivation of the expression for the 
deposited energy 

We start from a two-dimensional system at thermal 
equilibrium with no average current and we apply a gauge 
field of the form 

A(r,t) ^ f{t)e v \c+e^ r ' 2 + c_e^ r / 2 \ 2 e-^-^ 2 /^ 2 \ 

(CI) 

where e a is the unit vector along direction a. The deriv- 
able dimensionless envelope function /(£) is assumed to 
be zero for t < and to rapidly tend to zero for t — > +oo. 
The time-independent coefficients c± have the dimension 
of the square root of a momentum. 

We are interested in evaluating the energy change of 
the system from t — 0~ to t = +oo at the lowest order 
in c±. We work in Schrodingcr picture and we first use 
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the exact relations 
AE = 



~, at 



(C2) 



+°° r d 

dtj d 2 rA(r,t)--(j(r))(t), (C3) 



where the second equality comes from a time-dependent 
Hcllmann-Feynman theorem and a temporal integration 
by parts. Calculating (j(r))(f) by linear response theory 
gives 
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where x° x is the exact current susceptibility in real 
space taking into account the spatial inhomogeneity of 
the trapped cloud. Note that, contrarily to Eq. (|C2[) . 
Eas. (|U3l) . (|U4|) still hold when f(t) = for t < and 
has a discontinuous jump in t — 0. 

We now use the particular form (jCljl for A and con- 
sider the relevant limiting case qa S> 1, <7min(£, A) <C 1, 
where £ is the healing length of the gas and A is the 
thermal de Broglie wavelength. We also assume that a 
is much smaller than the radius of the trapped cloud, 
so that the density variation within a region of radius a 
around ro may be neglected. 

Within a local density approximation, we then re- 
place x ex with the susceptibility x of a spatially homo- 
geneous system with a density equal to the one of the 
trapped gas at position ro and with the same tempera- 
ture, x cx 0, r'; u) ~ x(* ~ r '; w )- 

This local density approximation leads to 
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where we have performed the change of variables r = 
R + u/2,r' = R u/2. As we work in the qa 1 
regime, we have for example 
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so that all the oscillating terms in R may be neglected. 
Introducing the Fourier transform of Xyyfcuj), which is 



an even function of k due to parity or rotational invari- 
ance, we obtain 
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The second contribution AE 2 comes from the spatially 
modulated gauge field at q, while the first contribution 
AEi is due to the non-modulated term which follows 
the broad Gaussian envelope. The expression (f24)) in 
the main text is obtained from (|C9|) by noting that the 
integration over k is effectively limited by the Gaussian 
factor to a small region of radius 1/er in which one is 
allowed to neglect the k-dependence of the susceptibility. 

Naively, one could guess that a necessary condition 
for the accuracy of our local density approximation is 
that the switch-off time of the gauge field 7 _1 is short 
as compared to the characteristic time 2R/v for the in- 
duced mechanical perturbation to cross the whole cloud, 
to be reflected by its boundaries and to turn back to the 
excitation zone where it can interfere with the excita- 
tion process, v being the fastest between the sound and 
thermal speeds in the cloud of radius R. 

This condition is actually sufficient, but not neces- 
sary within linear response theory. We now show for 
f(t) = e{t)e-^, as in Eq.flHD, that the 7 -> limit for 
the deposited energy scheme exists and coincides with the 
perturbation induced by the gauge field in the thermo- 
dynamic equilibrium state. As one can show by inserting 
the explicit form of the temporal Fourier transform of the 
gauge field into (|C4[) and performing the integral over ui, 
the deposited energy can be written in the form 



AE~^J d 2 rd 2 r' ^A Q (r,t = 



Ap(v',t = Q + )Re[xTp(ry;u = *7)] (CIO) 

where we have introduced the Kubo formula for the exact 
current-current susceptibility 
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in terms of the thermal equilibrium population nx of 
quantum state A. This quantity can be simply related 
to the susceptibility at thermodynamic equilibrium, 

X%{v,v') = lim ,Re^(r,r';, = i 7 )] + 

7r A (Ab a (r)|A')(A'b>(r')|A). (C12) 
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We recall that the thermodynamic susceptibility relates 
the mean current in a thermal equilibrium state at tem- 
perature T to the applied (weak) static gauge field via 
(j) = x th * A, where * is the spatial convolution prod- 
uct. In the present case of an interacting gas, one can 
safely expect that the second line in (|C12|) gives a neg- 
ligible contribution as there is no systematic degeneracy 
and the current operator j has no diagonal matrix ele- 
ments since the eigenstate wavefunctions may be taken 
real. Within linear response theory, we thus obtain 

AE -+ 1 f d 2 rd 2 /^ X t a h /3 (r,r'M Q (r,0 + )A /3 (r',0+). 

(C13) 

Since local density approximation in the thermodynamic 
equilibrium state is a standard procedure, we expect 
that it can be applied to evaluate the right-hand side 
of Eq. (|C13[) [73|]. As a consequence, within linear re- 
sponse theory, our deposited energy method (that mea- 
sures the left-hand side of Eq. (|C13p ) should be accurately 
described by the local density approximation down to the 
7 — > limit. 

Appendix D: Some results of linear response theory 
and the Bogoliubov expression of the deposited 
energy 

A system of time independent Hamiltonian H experi- 
ences, at times t > 0, a time dependent weak perturba- 
tion of Hamiltonian — ef(t)V, where e — > 0, the dimen- 
sionless time dependent factor f(t) is zero for t < and 
tends rapidly to zero for t — > +00, and the operator V is 
time independent. At time t = +00, the system is free 
again, with a mean energy modified by the perturbation. 
The question is to calculate the mean energy change to 
second order in e. 

Suppose first that, at t = 0~, the system is prepared in 
the eigenstate |A) of H of eigenenergy E x . The energy 
change SE between time and time +00 is 

5E = lim (i>(t)\(H Q - E x )\i>(t)) (Dl) 
t— 

where \ip(t)) is the system state vector at time t. The 
usual time dependent perturbation theory gives an ex- 
pansion in powers of e: 

\m) = \Mt)) + e\Mt)) + SIM*)) + ■■■ (D2) 
where \ipo(t)) = exp(— iE x t/h)\X), 

ft A T 

\M*)) = - / —e- lHo{t - T)/h f(T)Ve- lE * T/h \\), (D3) 
Jo ™ 

and the expression of higher order contributions is not 
needed. Using (H - E x )U' (t)} = and (V>o (t) | (H - 
E\) — 0, we find to second order in e: 



In this paper, f(t) = 0(i) cxp(— jt), with 7 > and 
8(t) is the Heaviside step function, see Eq. ([2T|) . Also, 
the system is prepared initially in a statistical mixture 
of eigenstates of Ho with a probability distribution tt\. 
After explicit integration of (|D3|) over r and then aver- 
age over I A) , the expression for the signal to be detected 
experimentally is 

Signal (V) = lim = 

e^O e 2 

The sum may be restricted to E\ 7^ Ey since the con- 
tributions with E\ = E\i are zero. This also shows that 
the signal has a finite limit for 7 — >• + . Note that in a 
thermal equilibrium state ir\ = Z^ 1 cxp(—E\/kBT), the 
signal is necessarily positive. 

The calculation of the noise on the experimental sig- 
nal can be performed along the same lines. One defines 
8E 2 = ma t -y +oo {ip(t)\(H - E x ) 2 \ip{t)) with the initial 
state vector 1-0(0)) = |A), and one finds after average 
over the initial state: 

[Noise (V)] 2 lim ^ = £ 7r A |(A|V|A')| 2 * 
^tt a [(A|V 2 |A)-(A|V|A) 2 ], (D6) 

A 

where the approximate equality is based on the assump- 
tion that there are no systematic degeneracies in the 
many-body spectrum. The e 2 scaling of the variance 
in Eq. (|D6p shows that the typical value of the energy 
change at the end of the excitation sequence is of order 
e. This scaling is to be contrasted with the e 2 one of the 
expectation value that is suggested by Eq. (|D5p . 

We now apply the general formula Eq. (|D5|) to the 
Bogoliubov analysis of subsection IIVBI In this case, 
e = £gau g c/2 and V = V q + V- q with 

V q = f d 2 re - Jy (r)=^^at +q a k (D7) 

where is the annihilation operator of a particle of the 
gas of wavevector k, and we have set q = qe x . In a 
translationally invariant system, the eigenstates |A) can 
be taken of well defined total momentum; as the action 
of V± q changes this total momentum by ±fiq, the two 
operators V q and V_ g cannot interfere in the signal and 
thus Signal (V) = 2 Signal (V q ). In terms of the annihi- 
lation operators of Bogoliubov quasiparticles, we can 
split 



5E~ lim e 2 (Mt)\(Ho-E x )\Mt))- (D4) 



(D8) 
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in terms of 



v f )= E ^(^k +q -m +q )4 + A( D9 ) 

k#0,-q 

k#0,-q 
' fife 

Mo-q r 



The primed sum indicates restriction of the sum over 
wavevectors such that k y > 0. In Bogoliubov theory, the 
cigenstatcs |A) may be taken in the form of Fock states 

of quasiparticles. Since changes the total number of 
quasiparticles by the amount n, the terms in the right 
hand side of Eq. (|D8|) cannot interfere in the signal and 



Signal (V,) = Signal (V^ 0) )+Signal(V( 2) )+Signal(V(- 2) ). 

(D12) 

(n) 

Thanks to the clever writing of Vq with the con- 
straint k y > [74| . there are no interferences in 
Eqs. (|D9ID10IDll|) between the different terms of the 
sums over k. As a result, the whole signal is the sum 
over the contribution of the different k's. The last trick 
is to express the ratios 7r\< / tt\ in terms of the mean oc- 
cupation numbers nk of the Bogoliubov modes of energy 
£k and to make use of the identity 



rt k + 1 
n k 



(D13) 



satisfied by the Bose law. A little bit of rewriting taking 
advantage of the relation Signal(V) = (N/m)f° s and of 
the remarks in [llj], finally leads to Eq. (|3Tj) . 

An alternative procedure is to calculate the current- 
current susceptibility (|C11[) for a spatially homogeneous 
system within the Bogoliubov theory, which for q _L e y 
gives in dimension d: 



Xra(q;w) 



f 



E 



7Jk - "k+q 



- £ k+q - 


e k — hiu 


-i0+ 
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- ek — hu 


-*0+ 


1 + 
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fq 


£k+q "+ 


£k — fojj 


-i0+ 


n k 
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£k — e k+q 



-hu- 



k#0,-q 

(^k+q^k - C/k+qVk+qC/kVk) 

(^k+q^k - U k+ qV k+ JJ k V k ) 

K +Cl v 2 - u k+q v k+tl u k v k ) 

U k +qV k+ qU k V k 
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From Eq. (|C4p one then recovers expression Eq. tpij) for 
the effective normal fraction /° ff . 

Another application of Eq. (|D5[) is to calculate the 
energy deposited by the external potential U (r, t) = 
(Uoe 1 ® ' r +c.c)®(t)e~ vt in the spatially homogeneous case, 



for -q > 0. This is useful for Appendix IA1 and Appendix [Bl 
to estimate the effect of undesired light-shifts. In second 
quantized form, and to leading order in Bogoliubov the- 
ory, one then has eV = N^Uq + V Q )po{b ] Q + 6_ Q ) + 
h.c.]. These terms do not interfere in Eq. (|D5p . For non- 
zero temperature, using (|D13|) . we then obtain a temper- 
ature independent result 



AE U ~ 2N\U \ 2 (U Q + V Q ) 2 Re 



f 



€ Q - ihrj 



(D15) 



Remarkably, this also allows to calculate the energy 
change AE g due to the switch-on-and-off of a spatially 
modulated coupling constant, Sg(r,t) = (<5<?oe iqr + 
c.c.)0(t)e~ 7t . For a spatially homogeneous system, to 
leading order of Bogoliubov theory, the relevant oper- 
ator is eV = NV 2 (U q + V q )[p5g<M + 6_ q ) + h.c], so 
that one can formally apply Eq. (|D15l) with Uq = pSgg. 
This can be applied to the variation of the coupling 
constant due to g aa =^ g a b in Eq. (|18p . In this case 
Sg = 2g[(g ab - g aa )/g aa }(n+n-*)(0+)/\n c \ 2 so that for 
q£ = 1 and j — > 0, 



16 mpg ( g ab - g aa 



AE„ 



AE 2 5/„ (hk c f V 9. 



(D16) 



For the values of Table Q] and \g a b — g aa \ < 0.1|p aa | as 
expected for 87 Rb, this gives AE g /AE 2 « 7 • 10~ 4 //„ 
which is negligible. 



Appendix E: The noise on the deposited energy 



In the section IIV C| while presenting the numerical 
results on the deposited energy measurement, we men- 
tioned that the statistical noise on the deposited energy 
was larger for smaller values of the gauge field amplitude 

Egauge- 

To understand this feature, it can be useful to rewrite 
the deposited energy for a single realization of the clas- 
sical field simulation in the form 

5E=- J d 2 rj(r,0) • A(r,0 + ) 

r+oo <• 

-J dt J d 2 r3(r,t)d t A(r,t). (El) 

The first term comes from the abrupt switch-on of the 
gauge field. For each realization, it is of order e ga ugc 
but averages to zero in the limit of an infinite number 
of realizations of the experiment as (j(r, 0)) = 0. In any 
actual calculation, an average over a finite number n rca i 
of realizations is taken, which gives a non-zero random 
value for 5E scaling as egauge/V^cai- 

The relevant signal (SE) is given by the second term 
in (jEip . obtained from the classical Hamiltonian iden- 
tity dH/dt = d t H. For small values of the gauge field 
switch-off rate 7, this term is of order 0(eg auge ) as in 
this limit j adiabatically follows the thermal equilibrium 
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value for the instantaneous value of the gauge field. As 
a result, the number of realizations that are needed to 
extract the signal out of the statistical noise due to the 
first term grows as |e gaU gc|~ 2 , which perfectly explains 
the numerical observation. 

Furthermore, it has been demonstrated by a number 
of recent cold atom experiments that noise is not al- 
ways just an hindrance but can be also a source of useful 
physical information [75l [76j . As a simple example, we 
consider here the amplitude of the noise on the energy 
that is deposited in the system at each realization of the 
experiment. This quantity is quantified by the average 
(SE 2 ) of the square of the deposited energy in the 7—^0 
limit. Looking at (|E1[) . it is immediate to see that in 
the small e ga uge limit the dominant contribution comes 
from the square of the first term, which suggests that the 
noise on the deposited energy is related to the variance 
of the instantaneous fluctuations of the current operator. 
From the fluctuation-dissipation theorem, this quantity 
can then be related to the normal fraction of the gas. 

This idea can be put on solid grounds by developing a 
full quantum calculation. The linear response theory cal- 
culation performed along the lines of Appendix [D] leads 
to the expression in Heiscnbcrg picture 



as usual as 



SVvH=/ dte lut {V{t)V(Q)). (E4) 



Under the assumption that most of the spectral weight of 
the V operator lies in the low-energy region hto <C ksT, 
we can approximate 1 — e ~ huj / k BT ^ ^/ftgj". This is a 
quite standard approximation of many-body theory and 
is generally accurate in the small q limit [78| . After a few 
manipulations, it leads to the general expression 

S vv (t = 0) = J — S vv (uj) -f^ 

= k B TRe[ X vv(u J = 0)}, (E5) 



where the equivalent of Eq.([26|) was used to obtain 
the last identity. An application of the fluctuation- 
dissipation relation (|E5|) to the susceptibility and the 
fluctuations of the mass current in liquid He can be found 
in [3. 

The link between the variance of the deposited energy 
and the normal fraction is immediately obtained by ap- 
plying (|E5[) to the specific operator 



lim <[#„(- 

7— s-0 



Ho(0)] 



([|d 2 rj(r).A(r,0+)] 2 ) (E2) 



where Hq is the unperturbed Hamiltonian, that is the 
Hamiltonian without coupling to the gauge field. This re- 
lation connects the variance of the quantum equivalent of 
the deposited energy to the instantaneous fluctuations of 
the current operator and confirms our expectation based 
on the classical field model. Moreover, combined with 
the fluctuation-dissipation theorem, it can be the start- 
ing point for another proposal to measure f n . 

For a generic hermitian operator V with vanishing 
diagonal matrix elements in the eigenbasis of Hq, the 
fluctuation-dissipation theorem of linear response theory 
relates the imaginary part of the susceptibility \ t° the 
Fourier transform of the correlation function Syy(£) = 
(V(t)V(0)), 



Im[xvvM] 
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2h 



1 



-hu)/k B T 



(E3) 



The Fourier transform Syv(^) of the correlation func- 
tion in the thermodynamic equilibrium state is defined 



V = - / cTrj(r) • A(r,0+) 



(E6) 



and isolating the contribution of the spatially modulated 
gauge field proportional to fl^tt^. In this way, using 
Eq. (|C10[) . one is led to the final expression 



KW 2 P_ ( Egauge x 2 

4 m 



fn 



1 



lim([flo(+oo) - -ffo(O)] 2 ), (E7) 



2k B T 7->o 

which demonstrates an alternative way of extracting the 
value of the normal fraction /„ from a measurement of 
the statistical variance of the deposited energy in a series 
of experiments. 

It is however crucial to note that a measurement of f n 
based on the relation (|E7|) requires taking expectation 
values of the Hamiltonian operator at different times. 
This may be experimentally challenging as it requires ei- 
ther a non-destructive measurement of the initial energy 
of the system at t = before switching on the gauge field, 
or a very precise a priori knowledge of its value in a sort 
of microcanonical ensemble 173 . 
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